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Abstract 



(— I , We develop a rigorous theory of non-local Hamiltonian structures, built on the no- 

tion of a non-local Poisson vertex algebra. As an application, we find conditions 
that guarantee applicability of the Lenard-Magri scheme of integrability to a pair of 
compatible non-local Hamiltonian structures. 
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Local Poisson brackets play a fundamental role in the theory of integrable systems. Recall 
that a local Poisson bracket is defined by (see e.g. [TF86]): 

(1.1) {ui{x),Uj{y)} = Hij{u{y),u'{y), . . .;dy)6{x - y) , 



where u = (ui, . . . , ui) is a vector valued function on a 1-dimensional manifold M, 5{x—y) 
is the (5-function: ^j^ff{y)S{x — y)dy = f{x), and H{d) = [Hij{d))._^^-^ is an ^ x ^ matrix 

^ . differential operator, whose coefficients are functions in u,u', . . . ,u^^\ One requires, in 

^ I addition, that (1.1) "satisfies the Lie algebra axioms". 

One of the ways to formulate this condition is as follows. Let V be an algebra of 

(n) 

differential polynomials in ui, . . . ,ue, i.e. the algebra of polynomials in u]: , i e I = 

{!,...,£}, n e Z+, with uf^^ = Ui and the derivation d, defined by du'f^ = u["''^^\ or its 
algebra of differential functions extension. The bracket (1.1) extends, by the Leibniz rule 
and bilinearity, to arbitrary f,geV: 

(1.2) {f(a:),g{y)}=Y, ^ f S f W ^^^."{^-(^)' (^)} " 

i,jeIm,neZ+ Su^ OUj 
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Applying integration by parts, we get the following bracket on V/5V: 

(1.3) {S/,b).J|.«(^)|. 

where J is the canonical quotient map V —>■ V/dV and |^ is the vector of variational 
derivatives ^ = YjneZ+ Y") • Then one requires that the bracket (1.3) satisfies 

the Lie algebra axioms. (The skewsymmetry of this bracket is equivalent to the skewad- 
jointness of H{d), but the Jacobi identity is a complicated system of non- linear PDE on 
its coefficients.) In this case the matrix differential operator H(d) is called a Hamiltonian 
structure. 

Given an element J/i e V/dV, called a Hamiltonian functional, the Hamiltonian equa- 
tion associated to H{d) is the following evolution equation: 

(1.4) ? = ^(5)?- 

^ ' dt ^ ' 6u 

For example, taking H{d) = d and h = ^{u^ + cuu"), we obtain the KdV equation: 
^ = 3uu' + cu'". 

Equation (1.4) is called integrable if is contained in an infinite dimensional abelian 
subalgebra A of the Lie algebra V/dV with bracket (1.3). Picking a basis {S^hn}neZ+ of 
A, we obtain a hierarchy of compatible integrable Hamiltonian equations: 

— = H{c)— , ne Z+ . 

dtn du 

An alternative approach, proposed in [BDSK09], is to apply the Fourier transform 
F{x,y) 1-^ ^j^dxe^^^~y^ F{x,y) to both sides of (1.2) to obtain the following "Master 
formula" [DSK06]: 

(1-5) {fx9}=Y. S T%(A + 5m.(A + 5)(-A-5)--^. 

i,jeIm,neZ+ (JUj OU^ 

For an arbitrary I x (. matrix differential operator H{d) this X-bracket is polynomial in 
A, i.e. it takes values in V[d], satisfies the left and right Leibniz rules: 

(1.6) {fxgh} = g{fxh} + h{fxg} , {fgxh} = {fx+sgUh + {fx+dh}^g , 

where the arrow means that X + d should be moved to the right, and the sesquilinerity 
axioms: 

(1.7) {dfxg] = -A{/a5} , {fxdg] = (A + d){fxg] . 

It is proved in [BDSK09] that the requirement that (1.3) satisfies the Lie algebra axioms 
is equivalent to the following two properties of (1.5): 

(1.8) {gxf} = -{f-x-eg} , 
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(1.9) {fx{9f.h}} = {g^Uxh}} + {{fx9]x+^.h} . 

A differential algebra V, endowed with a polynomial A-bracket, satisfying axioms (1-6)- 
(1.9), is called a Poisson vertex algebra (PVA). 

It was demonstrated in [BDSK09] that the PVA approach greatly simplifies the theory 
of integrable Hamiltonian PDE, based on local Poisson brackets. For example, equation 
(1.4) becomes, in terms of the A-bracket associated to H: 

and the Lie bracket (1.3) becomes 

It is the purpose of the present paper to demonstrate that the adequate (and in 
fact indispensable) tool for understanding non-local Poisson brackets is the theory of 
"non-local" PVA. 

We define a non-local X-bracket on the differential algebra V to take its values in 
V((A~^)), formal Laurent series in A~^ with coefficients in V, and to satisfy properties 
(1.6) and (1.7). The main example is the A-bracket given by the Master Formula (1.5), 
where H{d) is a matrix pseudodifferential operator. The only problem with this definition 
is the interpretation of the operator j^', this is defined by the geometric progression 

Property (1.8) of the A-bracket is interpreted in the same way, but the interpretation 
of property (1.9) is more subtle. Indeed, in general, we have {fx{9^h}} e V((A^^))((//~^)), 
but {g^Uxh}} e VHi^r^MX-')), and {{fx9}x+>^h} e V(((A + ^^)-i))((A-i)), so that all 
three terms of (1.9) lie in different spaces. Our key idea is to consider the space 

Vx,^ = V[[X-\|^^\{X + ^^)-']][\,^l], 

which is canonically embedded in all three of the above spaces. We say that a A-bracket 
is admissible if 

{f\{9fih}} e Va,^ , for aU f,g,heV. 

It is immediate to see that then the other two terms of (1.9) lie in Va.^ as well, hence 
(1.9) is an identity in Va,;^- 

We call the differential algebra V, endowed with a non-local A-bracket, a non-local 
PVA, if it satisfies (1.8), is admissible, and satisfies (1.9). 

For an arbitrary pseudodifferential operator H{d) the A-bracket (1.5) is not admis- 
sible, but it is admissible for any rational pseudodifferential operator, i.e. contained in 
the subalgebra generated by differential operators and their inverses. We show that, as 
in the local case (see [BDSK09]), this A-bracket satisfies conditions (1.8) and (1.9) if and 
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only if (1.8) holds for any pair Ui,Uj, and (1.9) holds for any triple Ui,Uj,Uk- Also, (1.8) 
is equivalent to skewadjointness of H{d). 

The simplest example of a non-local PVA corresponds to the skewadjoint operator 
H{d) = Then 

{uxu} = A"^ , 

and equation (1.9) trivially holds for the triple u,u,u. Note that (1.1) in this case reads: 
{u{x),u{y)} = d~^5{x — y), which is quite difficult to work with. 

The next example corresponds to Sokolov's skewadjoint operator [Sok84], H(d) = 
u'd~^ o u' . The corresponding A-bracket is 

{uxu] = u'—^u' . 
X + o 

The verification of (1.9) for the triple u,u,u is straightforward. 

A rational pseudodifferential operator H{d) is a Haniiltonian structure on V if the A- 
bracket (1.5) endows V with a structure of a non-local PVA (in other words H(d) should 
be skewadjoint and (1.9) should hold for any triple Ui,Uj,Uk)- 

Fix a "minimal fractional decomposition" H = AB^^. This means that A,B are 
differential operators over V, such that Kei A n Keri? = in any differential domain 
extension of V. It is shown in [CDSK12b] that such a decomposition always exists and 
that the above property is equivalent to the property that any common right factor of A 
and B is invertible over the field of fractions of V. Then the basic notions of the theory 
of integrable systems are defined as follows. A Hamiltonian functional (for H = AB~^) 

is an element J/i 6 V/dV such that = B{d)F for some F e V^. Then the element 
P = A{d)F is called an associated Hamiltonian vector field. Denote by J'{H) cz V/dV 
the subspace of all Hamiltonian functionals, and by %{H) a the subspace of all 
Hamiltonian vector fields (they are independent of the choice of the minimal fractional 
decomposition for H): 

(—\ ( Tm b\ c- Vmv 14{m= Ai B-^i Tin W 1/ 



Then it is easy to show that J-{H) is a Lie algebra with respect to the well-defined 
bracket (1.3), and ^{H) is a subalgebra of the Lie algebra with bracket [-P, Q] = 
DQ(d)P — Dp{d)Q, where Dp{d) is the Frechet derivative. 

A Hamiltonian equation, associated to the Hamiltonian structure H and a Hamilto- 
nian functional S^h e J^{H), with an associated Hamiltonian vector field P e ^{{H), is the 
following evolution equation: 

du „ 

Note that (1.10) coincides with (1.4) in the local case. The Hamiltonian equation (1.10) 
is called integrable if there exist linearly independent infinite sequences ^hn e J-{H) 
and Pn e H.{H), n e Z+, such that ^/iq = ^h, Pq = P, Pn is associated to ^hn, and 
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{^hm,S^hn} = 0, [Pm,Pn] = for all m,n e Z+. In this case we have a hierarchy of 
compatible integrable equations 



Having given rigorous definitions of the basic notions of the theory of Hamiltonian 
equations with non-local Hamiltonian structures, we proceed to establish some basic 
results of the theory. 

The first result is Theorem 5.1, which states that if H and K are compatible non-local 
Hamiltonian structures and K is invertible (as a pseudodifferential operator), then the 
sequence of rational pseudodifferential operators ifM = K, H^"'^ = {H o K~^)^~^ o H , n ^ 
1, is a compatible family of non-local Hamiltonian structures. (As usual [Mag78, MagSO] 
a collection of non-local Hamiltonian structures is called compatible if any their finite 
linear combination is again a non-local Hamiltonian structure.) This result was first 
stated in [MagSO] and its partial proof was given in [FF81] (of course, without having 
rigorous definitions). 

Next, we give a rigorous definition of a non-local symplectic structure and prove 
(the "well-known" fact) that, if S is invertible as pseudodifferential operator, then it is 
a non-local symplectic structure if and only if is a non-local Hamiltonian structure 
(Theorem 6.2). Since we completely described (local) symplectic structures in [BDSK09], 
this result provides a large collection of non-local Hamiltonian structures. We also estab- 
lish a connection between Dirac structures (see [Dor93] and [BDSK09]) with non-local 
Hamiltonian structures (Theorems 6.11 and 6.16). 

After that we discuss the Lenard-Margi scheme of integrability for a pair of com- 
patible non-local Hamiltonian structures, similar to that discussed in [Mag78, MagSO, 
Dor93, BDSK09] in the local case, and give sufficient conditions when this scheme works 
(Theorem 7.7 and Remark 7.S). 

Finally, we demonstrate on the example of the NLS hierarchy how this scheme works. 
The corresponding pair of compatible non-local Hamiltonian structures has been written 
down already in [MagSO]: 



In our next paper [DSK12] we go systematically over well known non-local bi-Hamiltonian 
structures and construct the corresponding integrable hierarchies of Hamiltonian PDE's. 

Throughout the paper, all vector spaces are considered over a field F of characteristic 
zero. 

We wish to thank Pavel Etingof and Andrea Maffei for useful discussions. 

2 Rational pseudodifferential operators 

2.1 The space Va,m 

Throughout the paper we shall use the following standard notation. Given a vector space 
V, we denote by V[X] the space of polynomials in A with coefficients in V, by V^[[A^^]] the 
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space of formal power series in A^^ with coefficients in V, and by V{{X^^)) = V[[X^^]][X] 
the space of formal Laurent series in A^^ with coefficients in V. 

We have the obvious identifications V[X, /i] = F[A][/x] = y[/x][A] and ^[[A^^, /i^^]] = 
= ^[[/^"^]][[A"^]]- However the space V{{X-^)){{fi-^) does not coincide 
with y((/i^^))((A^^)). Both spaces contain naturally the subspace 1/[[A~^, //~^]][A, 
In fact, this subspace is their intersection in the ambient space F[[A-^,^-^]] consisting 
of ah series of the form 0'm,n 

The most important for this paper will be the space 

^A,M:=^[[A-\M-\(A + A^)-i]][A,Ai], 

namely, the quotient of the F[A, /.i, i/]-module V\\X^^ , ^^"^ , i/^-'^]][A, fi, u] by the submodule 
(z^ — A — iJ,)V[[X^^ , fi^^ , u^^JllX, fi, u]. By definition, the space V^.^ consists of elements 
which can be written (NOT uniquely) in the form 

MNP 

(2.1) ^= S S S «™A-^"(A + /x)^ 

m=—co n=— oop=— 00 

for some M, N,P e 7j (in fact, we Cctn always choose ^ 0), cind. cirn,n,p ^ 

V. 

In the space y[[A~^, /i^^, i/^^]][A, /i, i^] we have a natural notion of degree, by let- 
ting deg(A) = deg(/i) = deg(z/) = 1. Every element A e V[[X^^ , fi^^ ,u^^]][X, fiju] 
decomposes as a sum A = (possibly infinite), where is a finite linear 

combination of monomials of degree d. Since — A — /i is homogenous (of degree 1), 
this induces a well-defined notion of degree on the quotient space Vx^^, and we denote by 
V^^^, for d e Z, the span of elements of degree d in Vx,^. If A e Vx^^ has the form (2.1), 
then it decomposes as ^ = Xidl-o)^^ > where A^*^^ e V^^ is given by 

A^''^= 2 a^,„,pA"^/i"(A + /i)^ 

(m+n+p=d) 

The coefficients arn,n,p £ ^ still not uniquely defined, but now the sum in A^'^^ is finite 
(since d — 2K ^ m^n^p ^ K := max(M, N, P)). Hence, we have the following equality 

^A'^ = nA±\/^±\(A + M)-Y, 

where, as before, the superscript d denotes the subspace consisting of polynomials in 
A-^, (A + /i)~^, of degree d. 

Lemma 2.1. The following is a basis of the space V^^ over V: 

X'^-'fi\ ieZ ; X'^+\X + i e Z>o = {1, 2, . . . } , 

in the sense that any element of the space V^^ can be written uniquely as a finite linear 
combination of the above elements with coefficients in V. 
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Proof. First, it suffices to prove the claim for d = 0. In this case, letting t = /i/A, 
the elements of ^ are rational functions in t with poles at and -1. But any such 
rational functions can be uniquely written, by partial fractions decomposition, as a linear 
combination of t*, with z e Z, and of (1 + with i e Z>o. □ 

Remark 2.2. One has natural embeddings of in all the spaces y((A~^))((/i~^)), 
y((/.-i))((A-i)), y((A-i))(((A + ^)-i)), y((/i-i))(((A + /.)-!)), F(((A + A.)-i))((A-i)), 
V{{{X + /x) ^)), defined by expanding one of the variables A, /x or = A + // in 

terms of the other two. For example, we have the embedding 

(2.2) i,,x: Vx,^^V{{X-'Mfi-')), 

obtained by expanding all negative powers of A + ;U in the region \fi\ > |A|: 

(2.3) v,a(a + ^^)—' = f; ("7 ') x'^^-'-' . 

Similarly in all other cases. Note that, even though V\^^ is naturally embedded in both 
spaces l^((A~^))((/x^^)) and y((/i^^))((A^^)), it is not contained in their intersection 
V[[X-\ f,-']][X, f,]. 

Lemma 2.3. // V is an algebra, then Va,^j is also an algebra, with the obvious product. 
Namely, if A{X, fi), B{X, n) e V\^^, then A{X, n)B{X, fj,) e Vx^^. More generally, if S,T : 
V are endomorphisms ofV (viewed as a vector space), then 

A{X + S,fi + T)B{X,fi)eVx,^, 

where we expand the negative powers of X + S and + T in non-negative powers of S 
and T, acting on the coefficients of B. 

Proof. We expand A and B as in (2.1): 

M N P 
A(A,/z)= 2 2 S «m,„,pA-M"(A + /i)^ 

■m=— CO n=— oc p=— oo 
M' N' P' 

B(A,/x)= 2 2 2 6^.,„.yA"^V"'(A + /i)^''. 

»n'=— cx) n'=— CO p'=— oo 

Using binomial expansion, we then get 

M+M' N+N' P+P' 

A{X + S,fi + T)BiX,^^)= 2 S S c^,n,p-A'^/u"(A + /i)^ 

m=— 00 n=— 00 j5=— oo 

where 

Cfh,n,p =2 2 S 

(m+m'— i=m) {n+n'—j—n) (p+p' —k=p) 



'rn' ,n' ,p' ) 



To conclude, we just observe that each sum in the RHS is finite, since, for example, in 
the first sum we have i = m + m' — fh, fh — M' ^ m ^ M and fh — M ^ m' ^ M' . □ 



7 



Lemma 2.4. Let V be a vector space and let U V be a subspace. Then we have: 
{AeVx,^\L^,xAeU{{X-^)){{fi-^))} 



.A,^AeC/((/x-i))((A-i))} 
ix,x+,AeU{i{X + ^^r')){{X-'))} 

i,^x+,AeU{{{x + ^l)-'M^l-^))} 



u 



\,fl ■ 



Proof. We only need to prove that {A e Vx,f, \ if,,xA e U{{X-^)){{fi-^))} cz C/^,^. Indeed, 
the opposite inclusion is obvious, and the argument for the other equalities is the same. 

Let A e V^^ be such that its expansion i^^xA e V{{X~^)){{iJ,~^)) has coefficients in 
U. We want to prove that A lies in Ux,^- By Lemma 2.1 A can be written uniquely as 



A 



N 



N 



ViX'^^'^-' 



-M 



with Vi,Wj e V. Its expansion in V{{X ^)){ip ^)) is 



N 



A 



2 ^^^'^ 



N 00 

EE 

i=i k=0 



-j-k 



Since, by assumption, tfj,^xA e U{{X ^)), we have 



Vi e U 



N 



-J 

-i 



-J 



Wj e U 



for 
for 



- M ^ i ^ 
^ i ^ iV, 



for i > N . 



From the first condition above we have that Vi lies in U for i < 0. Using the third 
condition, we want to deduce that wj lies in U for all 1 ^ j ^ N. It then follows, from 
the second condition, that Vi lies in U for i ^ as well, proving the claim. 

For i > and 1 ^ j ^ iV we have (rj.) = . Hence, we will be able to 

deduce that Wj lies in U for every j, once we prove that the following matrix 



P 



J-1 



has rank N. First, the sign (— 1)*"^-' does not change the rank of the above matrix. So it 
sufficies to prove that the matrices 



Th 



/i+l5Sis;/i+Af 



are non-degenerate for every h ^ 0. This is clear since the matrix Tq is upper triangular 
with ones on the diagonal, and, by the Tartaglia-Pascal triangle, and T/j+i have the 
same determinant. □ 
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2.2 Rational pseudodifferential operators 

For the rest of this section, let ^ be a differential algebra, i.e. a unital commutative 
associative algebra with a derivation 5, and assume that ^ is a domain. For a e A, we 
denote a' = d{a) and a^") = 5" (a), for a non negative integer n. We denote by /C the 
field of fractions of A. 

Recall that a pseudodifferential operator over A is an expression of the form 

N 

(2.4) A = A{d) = 2 fln^" ' (^n^A. 

n=— 00 

If Oat 7^ 0, one says that A has order N. Pseudodifferential operators form a unital 
associative algebra, denoted by A{{d^^)), with product o defined by letting 

(2.5) d''oa= 2 a^^^d""-^ , neZ,aeA. 

We will often omit o if no confusion may arise. 

Clearly, lC{{d~^)) is a skewfield, and it is the skewfield of fractions of A{{d^^)). If 
A e A{{d^^)) is a non-zero pseudodifferential operator of order as in (2.4), its inverse 
A^^ e lC{{d^^)) is computed as follows. We write 

-1 

A = aN{l+ 2 «^V+7v5")5^, 

n=— 00 

and expanding by geometric progression, we get 

00—1 ^ 

(2.6) ^-1 = o 2 ( - 2 a^ia„+Ar5") o a],' , 

fc=0 n=-oo 

which is well defined as a pseudodifferential operator in /C((5^^)), since, by formula (2.5), 
the powers of d are bounded above by —N, and the coefficient of each power of 5 is a 
finite sum. 

The symbol of the pseudodifferential operator A(d) in (2.4) is the formal Laurent series 
= 1j^=-oo^"'^" ^ "^(('^ "'^))' where A is an indeterminate commuting with A. We 
thus get a bijective map A{{d^^)) —>■ ^((A^^)) (which is not an algebra homomorphism) . 
A closed formula for the associative product in A{{d~^)) in terms of the corresponding 
symbols is the following: 

(2.7) {AoB){X) = A{X + d)B{X). 

Here and further on, we always expand an expression as (A + 5)", n e Z, in non-negative 
powers of d: 
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Therefore, the RHS of (2.7) means S^Iq (7)«m^^r A^+^^J'. 

The algebra A{{d^^)) contains the algebra of differential operators A[d] as a subal- 
gebra. 

Definition 2.5. The field lC{d) of rational pseudodifferential operators is the smallest 
subskewfield of K,{{d-^)) containing A[d]. We denote A{d) = IC{d) n A{{d-'^)), the 
subalgebra of rational pseudodifferential operators with coefficients in A. 

The following Proposition (see [CDSK12, Prop. 3. 4]) describes explicitly the skewfield 
K.{d) of rational pseudodifferential operators. 

Proposition 2.6. Let A he a differential domain, and let K, he its field of fractions. 

(a) Every rational pseudodifferential operator L e K,{d) can he written as a right (resp. 
left) fraction L = AS~^ (resp. L = S^^A) for some A,Se A[d] with S 0. 

(h) Let L = AS^^ (resp. L = S^^A), with A,Se A[d], he a decomposition of L e JC{d) 
such that S has minimal possihle order. Then any other decomposition L = AiS^^ 
(resp. L = S^^Ai), with Ai,Si e A[d], we have Ai = AK , Si = SK (resp. 
Ai = KA, Ai = KS), for some K e /C[5]. 

2.3 Rational matrix pseudodifferential operators 

Definition 2.7. A matrix pseudodifferential operator A e Mat^x^ ^(("5^^)) is called ra- 
tional with coefficients in A if its entries are rational pseudodifferential operators with 
coefficients in A. In other words, the algebra of rational matrix pseudodifferential oper- 
ators with coefficients in A is Mat^x^ -4.(5). 

Let M = (^AijB^^^ . be a rational matrix pseudodifferential operator with coeffi- 
cients in A, with Aij,Bij e A[d]. By the Ore condition (see e.g. [CDSK12]), we can find 
a common right multiple B e A[d] of all operators Bij, i.e. for every i,j we can factor 
B = BijCij for some Cij e A[d]. Hence, AijB~^ = AijB~^, where Aij = AijCij. Then, 
the matrix M can be represented as a ratio of two matrices: M = A{B1)^^. Hence, 

MatexeAid) = f^A{Biy 

However, in general this is not a representation of the rational matrix M in "minimal 
terms" (see Definition 2.11 below). 

We recall now some linear algebra over the skewfield /C((5~^)) and, in particular, 
the notion of the Dieudonne determinant (see [Art57] for an overview over an arbitrary 
skewfield). 

An elementary row operation of an £ x £ matrix pseudodifferential operator A is either 
a permutation of two rows of it, or the operation T{i,j',P), where 1 ^ i j ^ m and 
P e /C((5~^)), which replaces the j-th row by itself minus i-th. row multiplied on the left by 
P. Using the usual Gauss elimination, we can get the (well known) analogues of standard 
linear algebra theorems for matrix pseudodifferential operators. In particular, any matrix 



A e Mat i^iA[d], B e A[d], \ 
A,jB-^eA{{d-^)) \/i,j J • 
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pseudodifferential operator A e Matmxe can be brought by elementary row 

operations to a row echelon form. 

The Dieudonne determinant of a ^ e Mat^x^ has the form det^ = c^*^, 

where c e ^ is an indeterminate, and (i e Z. It is defined by the following properties: 
det A changes sign if we permute two rows of A, and it is unchanged under any elemen- 
tary row operation T{i,j;P) defined above, for aribtrary i ¥= j and a pseudodifferential 
operator P e }C({d~^)); moreover, if A is upper triangular, with diagonal entries An of 
order and leading coefficoent aj, then 



It was proved in [Dic43] (for any skewfield) that the Dieudonne determinant is well 
defined and det{AoB) = (det A) (det B) for every ixi matrix pseudodifferential operators 
^,BeMat,xf/C((5-i)). 

The Dieudonne determinant gives a way to characterize invertible matrix pseudodif- 
ferential operators, thanks to the following well known fact (see e.g. [CDSK12, Prop.4.3]): 

Proposition 2.8. LetD be a subskewfield of the skewfield K{{d^^)) , and let A eMat^x^^- 
Then A is invertible in Mat^x^^ if o-nd only if Aet A ^ 0. 

Corollary 2.9. Let A e Mat£x£ ^^((5^"'^)) be a matrix with deiA 7^ 0. Then A is a 
rational matrix if and only if A^^ is. 

Proof. It is a special case of Proposition 2.8 when V is the subskewfield /C(5) JC{{d^^)) 
of rational pseudodifferential operators. □ 

Remark 2.10. It is proved in [CDSK12] that, if A e Mati^i A{{d'^)) then we have 
det A = G^'^, with c e A, where A is the integral closure of A. Furthermore, if c is an 
invertible element of A, then the inverse matrix A''^ lies in Mat^xi A{{d^^)). 

Definition 2.11 (see [CDSK12b]). Let H e Mat^x^^(5) be a rational matrix pseudodif- 
ferential operator with coefficients in the differential field /C. A fractional decomposition 
H = AB^^, with A,Be Matixi}C[d], det B ^ 0, is called minimal if deg^deti? is 
minimal (recall that it is a non- negative integer). 

Proposition 2.12 ([CDSK12b]). (a) A fractional decomposition H=AB^^ of a rational 
matrix pseudodifferential operator H e Mat£x£ A^(5) is minimal if and only if 



in any differential field extension o//C. 

(b) The minimal fractional decomposition of H exists and is unique up to multiplication 
of A and B on the right by a matrix differential operator D which is invertible in the 
algebra Mat^xf 




(2.9) 



Ker A n Ker S = 
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Remark 2.13. In the case i = 1 the fractional decomposition H = AB^^ e IC{d), is 
minimal if and only if the the differential operators A,B e IC[d] have no right common 
divisor of order greater than (i.e. the right greatest common divisor of A and B is 1). 

Remark 2.14. Let ^ be a differential domain, and let /C be its field of fractions. A 
fractional decomposition H = AB~^ of H e Matixe ^[S] over /C can be turned into a 
fractional decomposition over A by clearing the denominators of A and B. Hence, a 
minimal fracitonal decomposition over A is also minimal over /C. 

3 Non-local Poisson vertex algebras 

3.1 Non-local A-brackets and non-local Lie conformal algebras 

Let i? be a module over the algebra of polynomials F[5]. 

Definition 3.1. A non-local X-bracket on is a linear map {• a •} : R ^ R{{X~^)) 
satisfying the following sesquilinearity conditions (a, 5 e R): 

(3.1) {daxb} = -X{axb} , {axdb} = (A + d){axb} . 

The non-local A-bracket {-a"} is said to be skewsymmetric (respectively symmetric) if 
(a, b e R) 

(3.2) {bxa} = -{a_x-eb} ( resp. = {a^x^gb}^ . 

The RHS of the skewsymmetry condition should be interpreted as follows: if {axb} = 
S^^^ooCnA", then 

N Af 00 / X 

{a^x-sb}= 2 (-A-5)"c„= 2 S (^)(-l)'^(5'c„)A'^-' 

n=— 00 n=— oofc=0 ^ ^ 

N N-m / : l.\ 

= S ( S r, )(-i)™^'^(5W))a™. 

m=-co fc=0 ^ ^ 

In other words, we move —A — 5 to the left and we expand in non negative powers of d 
as in (2.8). 

In general we have {ax{b^c}} e R{{X^^)){{fi^^)) for an arbitrary A-bracket {-a'}- 
Recall from Section 2.1 that Rx,^i can be considered as a subspace of R{{X^^)){{fi^^)) 
via the embedding l^^x- 

Definition 3.2. The non-local A-bracket {-a •} on R is called admissible if 

{axib^c}} e Rx,^ 'ia,b,ceR. 

Remark 3.3. If {-a •} is a skewsymmetric admissible A-bracket on R, then {b^{axc}} e 
Rx,ii and {{axb}x+ij.c} e Rx,ij. for all a,b,c e R. Indeed, the first claim is obvious since 
Rx,fi = Rfi,x- For the second claim, by skewsymmetry {{axb}x+fic} = -{c^x-fj.-d{axb}}, 
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and by the admissibility assumption {ci,{axb}} e Rx^u- To conclude it suffices to note that 
when replacing by —A — /i — 5 in an element of Rx^^ = R[[X^^, u^^, (A + !/)^^]][A, u], we 
have that is expanded in negative powers of X+jj, and (A+i^)^^ is expanded in negative 
powers of fi. As a result, we get an element of R[[X~^ , , (A + /i)^^]][A, /i] = Rx,^- 

Definition 3.4. A non-local Lie conformal algebra is an F[5]-module R endowed with 
an admissible skewsymmetric A-bracket {• a ■} ^ R(S) R^ R{{X~^)) satisfying the Jacobi 
identity (in Rx,fj.)- 

(3.3) {ax{b^c}} - {b^{axc}} = {{axb}x+^,c} for every a,b,ce R. 

Example 3.5. Let R = {¥[d] ® F) © FC, where V is a vector space with a symmetric 
bilinear form (• | •). Define the (non-local) A-bracket on R by letting C be a central 
element, defining 

{axb} = {a\b)CX-^ foTa,beV, 

and extending it to a A-bracket on R by sesquilinearity. Skewsymmetry for this A-bracket 
holds since, by assumption, (• | •) is symmetric. Moreover, since any triple A-bracket is 
zero, the A-bracket is obviously admissible and it satisfies the Jacobi identity. Hence, we 
have a non-local Lie conformal algebra. 

3.2 Non-local Poisson vertex algebras 

Let V be a differential algebra, i.e. a unital commutative associative algebra with a 
derivation d : V ^ V. As before, we assume that V is a domain and denote by K, its field 
of fractions. 

Definition 3.6. (a) A non-local X-bracket on the differential algebra V is a linear map 
{•A'} : V ® V -* V((A^^)) satisfying the sesquilinearity conditions (3.1) and the 
following left and right Leibniz rules: 

{axbc} = b{axc} + c{axb} , 

{abxc} = {ax+8c}^b+ {bx+dc}^a. 

Here and further an expression {ax+db}^c is interpreted as follows: if {axb} = 
Xi^=-oo then {ax+db}^c = Cn(A + d)^c, where we expand (A + d)^c 

in non-negative powers of d as in (2.8). 

(b) The conditions of (skew) symmetry, admissibility and Jacobi identity for a non-local 
A-bracket {■ x •} on V are the same as in Definitions 3.1, 3.2 and 3.4 respectively. 

(c) A non-local Poisson vertex algebra is a differential algebra V endowed with a non- 
local Poisson X-bracket, i.e. a skewsymetric admissible non-local A-bracket, satisfying 
the Jacobi identity. 

Example 3.7 (cf. Example 3.5). Let V = \ i = 1, . . . ,£,n e Z+] be the algebra of 

diffenertial polynoamials in £ differential variables ui, . . . ,U£. Let C = [Cij)- be an 
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^ X £ symmetric matrix with coefficients in F. The following formula defines a structure 
of a non-local Poisson vertex algebra on V: 



For example, {ui\Uj} = CijX^^ but, for P,Q ^ . . . , U(\ cz V, we get an infinite formal 
Laurent series in A^^: 



00 



vSS(-l)"||(^"f)^-"-'^V((A-)). 

jj = ln=0 J * 

We will prove that this is indeed a non-local Poisson A-bracket in the next section, where 
we will discuss a general construction of non-local Poisson vertex algebras, which will 
include this example as a special case (see Theorem 4.1). 

Proposition 3.8. Let {■ x ■} be a non-local Poisson vertex algebra structure on the dif- 
ferential domain V. Then there is a unique way to extend it to a non-local Poisson vertex 
algebra structure on the differential field of fractions K, and it can be computed using the 
following formulas (a, b e /C\{0} ): 

(3.5) {axb-'} = ^b-^{axb} , {a^^b} = ~{ax+ebUa'^ . 

Proof It is straightforward to check that formulas (3.5) define a non-local A-bracket on 
the field of fraction /C, satysfying all the axioms of non-local Poisson vertex algebra. In 
particular, admissibility of the A-bracket can be derived from Lemma 2.3. The details of 
the proof are left to the reader. □ 

Thanks to Proposition 3.8 we can extend, uniquely, a non-local Poisson vertex algebra 
A-bracket on V to its field of fractions /C. The following results are useful to prove 
admissibility of a non-local A-bracket. 

Lemma 3.9. Let V be a differential algebra, endowed with a non-local X-bracket {■ x }• 
Assume that V is a domain, and let JC be its field of fractions. Let S = e 
Mat£x£ (^(((^ "'^))) be an invertible £ x £ matrix pseudodifferential operator with coeffi- 
cients in /C. Letting Sij = Xi^=-oc ^ij;nd^, the following identities hold for every a e K, 
and i, j e L: 

I N 

.3_g^ {ax{S-%{^i)} = -22 V,A(5-^).r(A + 11 + 3) 

{axsrt;n}{f^ + 5)" (^'^ )<, (/.) e /C( (A"! )) ((//-^ )) , 
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and 



I N 

((A + dT{S-%,{\)) .A,A+,.(5*-^)„(/.) e /C(((A + m)-^))((A-1)) , 

where l^ x and iA,A+/i ^'"e as in (2.3). /n equation (3.7) 5* denotes the adjoint of the 
matrix differential operator S (its inverse being (S^^)*). 

Proof. The identity 5 o = 1 becomes, in terms of symbols, 



t=i 



Taking A-bracket with a, we have, by sesquihnearity and the (left) Leibniz rule. 



e 

i N 
f=l n=— CO 

e N 

= 2 2 {axSrt;nKl^ + dnS-%{^) 
t=l 00 

e 

+ 2 V,A'5rt(A + /u + 5){aA(S-^)tj(^t)} . 

Note that if_i^xS{X + fi + d) is invertible in Mat^x^ (^[^^]((A^^))((/^i^^))), its inverse being 
i'fi,xS^^{X + fi + d). We then apply L^,x{S^^)ir{X + /i + 5) on the left to both sides of the 
above equation and we sum over r = to get 

e 

i e N 

= -SS S >^^^AS-')^AX + |^ + ^){axSrt■nK^^ + ^nS-%{p), 

r=l t=l n=~oo 



proving equation (3.6). 



15 



Similarly, for the second equation we have, by the right Leibniz rule, 

(. 



= 2{^rt(A + a)(5-i)t,(A)A+^a} 

i N 

= 2 2 {srtA>^ + dnS~%{X)x+,a} 

e N 

e N 

+ 2 2 {iS~^)tjWx+f,+da}^iX,X+^L{X - X - fl - d)'^Srt;n 
t=l n=— 00 

e N 

= 2 2 {'rt;n,^,+e^}^{X + dnS-%{X) 
t=l n=— 00 

e 

+ 2 {(^^^)ijWx+fi+da}^L.X,X+fiS*r{fJ') ■ 



We next replace in the above equation fi (placed at the right) hy fi + d, and we apply the 
resulting differential operator to Lx^\+^{S*^^)ri{fJ-)- As a result we get, after summing 
over r = 1, . . . 

I i N 

= -222 {^^•*;-A+M+a4^((A + a)"(5-i)t,(A)).A,A+M(5*"')«(^) , 

r=l i=l n=0 

proving equation (3.7). □ 

Corollary 3.10. Let V be a differential algebra, endowed with a non-local X-bracket 
{■ X •}. Assume that V is a domain, and let JC be its field of fractions. Let S = e 

Mat^x^ (^[^]) have non-zero Dieudonne determinant. Then the following identities hold 
for every a e IC and i, j e /.■ 



{a,{S-%{f,)} 

(3.8) =-22 (^"')-(^ + + S){axSrt;nKl^ + dnS-%{fl) E IC 



and 



(3.9) 



{iS-%{X)x+,a} 



e N 



2 Y^{srt-,n,+^+e^U{i^ + SriS-%W){S'-%.{fi) e/CA,,. 



r,t=l n=0 



where Sij = Xi^=o^u;"^"- 



16 



Proof. It is immediate from equations (3.6) and (3.7). 



□ 



Corollary 3.11. Let V be a differential algebra, endowed with a non-local X-bracket 
{■\}- Assume that V is a domain, and let JC he its field of fractions. Let A e V{d) = 
/C(5) n V{{d^'^)) be a rational pseudodifferential operator with coefficients in V. Then 
{axA{ii)} and {A{X)x+^a} lie in Vx^fj, for every a e V. 

Proof. First, note that if the pseudodifferential operators A,Be IC{{d^^)) satisfy the 
conditions 



and both terms in the RHS lie in JCx^fj. by the assumption on A and B, thanks to Lemma 
2.3. Similarly, by the right Leibniz rule. 



and both terms in the RHS lie in /C^,/^ (rather in the image of /C^.^ in ^(((-^+A^) "'^))((-^^^)) 
via iA,A+/^) by Lemma 2.3. By Corollary 3.10 we have that, if 5 e V[d], then {axS^^{fi)} 
and {S~^{X)x+^a} lie in /Ca,^ for all a e K. Hence, by Definition 2.5 and the above 
observations, we get that, if ^ e V{d) = K,{d) n V((5~^)), then {axA{^)} and {A{X)x+^l(''} 
lie in fCx,^ for all a e fC. On the other hand, if a e V, we clearly have {axA{^)} e 
V((A-i))((/x-i)) and {A{X)x+^,a} e V(((A+/x)-i))((A-i)). The claim follows from Lemma 
2.4 apphed to V = JC and U = V. □ 

Remark 3.12. In the case when S e V(5) is a rational pseudodifferential operator with 
coefficients in V, thanks to Corollary 3.11, we can drop l^ x and iA,A+^ respectively from 
equations (3.6) and (3.7), which hold in the space Vx,^- 

4 Non-local Hamiltonian structures 
4.1 Algebras of differential functions 

Let Re = F[m-"^ \ i e I,n e Z+] be the algebra of differential polynomials in the i 

variables Uj, i e / = {!,...,£}, with the derivation d defined by d{u\^'') = u^^'^^K The 
partial derivatives — are commuting derivations of Rg, and they satisfy the following 

commutation relations with d: 




{{A o 5)(A)A+^a} = {^(A + d)B{X)x+,,a} 

= {B{X)x+^,+pa}^Lx,x+^,A*{^^) + {A{X + ^)x+^,+^a}^B{X) 



(4.1) 



d 



d 



d 



(the RHS is if n 
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Recall from [BDSK09] that an algebra of differential functions is a differential algebra 
extension V oi Rg, endowed with commuting derivations 

extending the usual partial derivatives on Ri, such that only a finite number of — are 

non-zero for each / e V, and such that the commutation rules (4.1) hold on V. Examples 
other than Ri itself are any localization of Rg by a multiplicative subset or any algebraic 
extension of Rg. 

For / E V, as usual we denote by J/ the image of / in the quotient space V/dV. Recall 
that, by (4.1), we have a well-defined variational derivative ^ : V/dV V®^, given by 

neZ+ du] 

Note that if the algebra of differential functions V is a domain, then its field of 
fractions JC is again an algebra of differential functions in the same variables ui, . . . ,U£, 
with the maps — 7^77 : /C ^ /C defined in the obvious way. When V = /C we call it a field 

of differential functions. 

We denote by C = {c e V | 5c = 0} c V the subalgebra of constants, and by 



for all i e /, n e Z+ [ c V 



the subalgebra of quasiconstants. It is easy to see that C ^ T . Given / e V which is not 
a quasiconstant, we say that is has differential order N if —7777 for some i e I, and 

-^{w) = every j e I and n > N. We also set the differential order of a quasiconstant 

element equal to —00. We let Vat be the subalgebra of elements of differential order at 
most N. This gives an increasing sequence of subalgebras 

C^T = V_oo c Vo c Vi c • • • cz V , 

such that BVn Vn+i- 



4.2 Construction of non-local Hamiltonian structures 

Let V be an algebra of differential functions in the variables ui, . . . ,Ui. Assume that V is 
a domain, and let /C be the corresponding field of differential functions of fractions. Let 
H = {Hij^^.^j e Mat£x£ V((5^^)) be an £ x ^ matrix pseudodifferential operator over V, 
namely 

TV 

Hi,= 2 ^u;n5" eV((5-^)), i,jel. 

n=— 00 
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We associate to this matrix H a non-local A-bracket on V given by the following Master 
Formula (cf. [DSK06]) 

(4.2) {fx9}H= 2 -^{X + drH,,{X + d){-X-dr^ eV{{X-')). 

i,jel Su) du\ 

In particular 

(4.3) {uixUj}H = Hji{X) , i,j e I . 

The following result gives a way to check if a matrix pseudodifferential operator 
H e MatixiV{{d~^)) defines a structure of non-local Poisson vertex algebra on V. The 
analogous statement in the local case was proved in [BDSK09]. 

Theorem 4.1. Let V be an algebra of differential functions, which is a domain, ane let 
K, be its field of fractions. Let H e Mat^ V((c'^"'^)). Then: 

(a) Formula (4.2) gives a well-defined non-local X-bracket on V. 

(b) This non-local X-bracket is skewsymmetric if and only if IL is a skew-adjoint matrix 
pseudodifferential operator. 

(c) If H = (^Hij^^.^j e Mat^x^ V((5) is a rational matrix pseudodifferential operator with 
coefficients in V, then the corresponding non-local X-bracket {■ x'}h : V x V — > 
V((A~^)) (given by equation (4.2)j is admissible. 

(d) Let H = (^Hij^._.^j e Mat^x^ V(c/) be a skewadjoint rational matrix pseudodifferential 
operator with coefficients in V. Then the non-local X-bracket {■ x ■}h defined by (4.2) 
is a Poisson non-local X-bracket, i.e. it satisfies the Jacobi identity (3.3), if and only 
if the Jacobi identity holds on generators (i,j,k e I): 

(4.4) {uix{Uj^Uk}H}H - {uj^{uixUk}H}H - {{UixUj}Hx+^,Uk}H = , 

where the equality holds in the space Vx,f_i- 

Proof. For the proofs of (a), (b) and (d) one does the same computations as in the proof 
of [BDSK09, Thm.1.15] for the local case. So, we only prove part (c). Let a, f,g e V. By 
the Master Formula (4.2) and the left Leibniz rule, we have 

{ax{f,9}H}H= 2 {ax^{f. + drH,,{t, + d){-f,-dr^}^ 
i,jel du) ou\ 

i,jel OUj OUi 



+ S -4^^{\ + ^^ + ^T{axH,,{^i + ^)]H^{-^^-^Y 

ijel ou\ du) 



m,nsZ+ 

2 -% (A + + dTH,,{X + ^ + d){-X-^- dr{ax^}^ . 
du); du) 



m,neZ 
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All sums in the above equations are finite. Therefore, all three terms in the RHS lie in 
Va.^, thanks to Corollary 3.11 and Lemma 2.3. □ 



Definition 4.2. Let V be an algebra of differential functions, which is a domain. A 
non-local Hamiltonian structure on V is a skewadjoint rational matrix pseudodifferential 
operator with coefficients inV, H = (^Hij)..^^ e Matfx£V(5), satisfying equation (4.4) 
for every i, j, /c e /. 

Remark 4.3. It is easy to show that, if L e KL{d) is a rational pseudodifferential operator, 
then it can be expanded as 

00 N 

(4.5) ^=2 2 S Pid-^°P2d-^o...d-^op,d^, 

s=l n=Opi,...,psSVM 
(finite sum) 

for some fixed M,N e Z+. To see this, write L = AS^^, where A, S e V[d] and 
S = 2^=0 ^n^" non-zero leading coefficient sn, and expand using geometric 
progression: 

00 

(4.6) = 2 ( - ^'N^N-id-' s-^hod-"") s],' . 

On the other hand, it is not hard to see that if L admits an expansion as in (4.5), 
then {axL{fi)}H e -f^A,^ for every a e IC and every matrix pseudodifferential operator H. 
As a consequence, if all the entries of a matrix pseudodifferential operator H admit an 
expansion as in (4.5), then the corresponding A-bracket {■ \ •}h on /C is admissible. 

Remark 4.4. It is claimed in the literature (without a proof) [DN89] that, in order to 
show that a skewadjoint operator H defines a (local) Hamiltonian structure, it suffices to 
check the Jacobi identity for the Lie bracket {• , ■}h = {• a •}_ff|;^^o V/c^V on triples of 
elements of the form ^fui, where / e is a quasiconstant. This is indeed true, provided 
that the algebra of quasiconstants J- is "big enough", by the following argument. By a 
straightforward computation, using the Master Formula, we get 

ilfUi, ilgUj, '\hUk}H}H - {'\gUj,{lfUi, '\hUk}H}H 
-{{yUi,S^gUj}H,^hUk}H = ^h(^{uix{Uj^Uk}H}H 

-{uj^{uixUk}H}H - {{Ui\Uj}Hx+f^Uk}H'^{\x=df){\fi=0g) ■ 

Clearly, this is zero for all f,g,heT and all i,j, A; e / if and only if if is a Hamiltonian 
structure, provided that the algebra T satisfies the following non-degeneracy conditions: 

(i) if ^ha = for some a e V and all h e J^, then a = 0, 

(ii) if P{d)f = for some differential operator P e V[d] and for all f e F, then P = 0. 

Obviously, J- fulfills these conditions if it contains the algebra of polynomials ¥\x\. Often 
in the literature this criterion is used also for non-local Hamiltonian structures, which 
does not seem to have much sense, since in the non-local case V/dV does not have a Lie 
algebra structure. 
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4.3 Examples 

Example 4.5. Let V be any algebra of differential functions in i differential variables, 
with subalgebra of quasiconstants J-" c V. Any skewadjoint rational matrix pseudod- 
ifferential operator with quasiconstant coefficients, H = (^Hij(d)^ ..^j e Mat^-^i T(d), is 
a Hamiltonian structure. Indeed, by askewadjointness of H the A-bracket {-a"}// is 
skewsymmetric, and by the Master Formula, all triple A-brackets are zero. Note that, if 
H e MatixeJ'{{d^^)) is skewadjoint, even if it is not a rational matrix, the corresponding 
A-bracket {■ \-}h is still admissible, hence it defines a non-local Poisson vertex algebra 
on V. 

In the special case when Hij{\) = CijX^^, and C = {cij)jj^i is a symmetric matrix 
with constant coefficients, we recover the non-local Poisson vertex algebras from Example 
3.7. When C if a symmetrized Cartan matrix or extended Cartan matrix of a simple Lie 
algebra, we get the Hamiltonian structure for a Toda lattice (see [Fr98]). 

Example 4.6. The following three operators form a compatible family of non-local Hamil- 
tonian structures (i.e. any their linear combination is a non-local Hamiltonian structure) 
on the algebra Ri = ¥[u, u' , u" , . . .] of differential polynomials in one variable: 

(i) Ki = d (GFZ Hamiltonina structure), 

(ii) K^i = d^^ (Toda non-local Hamiltonian structure), 

(iii) H = u'd^^ o u' (Sokolov non-local Hamitonian structure), 

First, any linear combination over C of Ki and K^i is a non-local Hamiltonian structure, 
as discussed in Example 4.5. Next, it is easy to show (cf. [BDSK09, Example 3.14]) that 
is a symplectic structure on the field of fractions /Ci = Fraci^i, known as the Sokolov 
symplectic structure, [Sok84]. Hence, by Theorem 6.2 below, we deduce that H is a non- 
local Hamiltonian structure. To conclude that Ki, K-i, H form a compatible family, it 
suffices to check that 

where H{\) = u'{d + X)^^u' e V((A^^)). This is straightforward. 

Example 4.7. Dorfman non-local Hamiltonian structure on the algebra of differential 
polynomials Ri = F[u, u', u", . • • ] is: 

H = d-^o u'd-^ o u'd-^ . 

One easily shows (cf. [BDSK09, Example 3.14]) that is a symplectic structure 

on the field of fractions /Ci = Fraci?i, known as Dorfman symplectic structre, [Dor93], 
hence H is indeed a non-local Hamiltonian structure. Furthermore, one can show, by a 
lengthy calculation, that Sokolov's and Dorfman's non-local Hamiltonian structures are 
compatible. 

Example 4.8 (cf. [Dor93]). Another triple of compatible non-local Hamiltonian structures 
on Ri = F[n, u' ,u", . . .] is: 
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(i) Ki = d (GFZ Hamiltonian structure), 

(ii) K_i = d^^ (Toda non-local Hamiltonian structure), 

(iii) H = d^^ o u' + u'd^^ (potential Virasoro-Magri non-local Hamiltonian structure). 

Example 4.9 (cf. [MagSO]). There is yet another triple of compatible non-local Hamilto- 
nian structures on Ri = F[n, n', u" , . . . ]: 

(i) Ki = d (GFZ Hamiltonian structure), 

(ii) = d\ 

(iii) H = do ud^^ o ud (modified Virasoro-Magri non-local Hamiltonian structure). 

Example 4.10 (cf. [Mag 78, MagSO]). The following is a triple of compatible non-local 
Hamiltonian structures on R2 = F[n, v, u' , v' , . . .]: 

(i) Ki = dl (GFZ Hamiltonian structure), 

-1 



(iii) H = ( -_i I (NLS non-local Hamiltonian structure). 

' —UO O V UO O U ' 

5 Constructing families of compatible non-local 
Hamiltonian structures 



As in the previous sections, let V be an algebra of differential functions in the variables 
ui, . . . , ue, we assume that V is a domain, and we let IC be its field of fractions. As in the 
local case, two non-local Poisson vertex algebra A-brackets on V (respectively two non- 
local Hamiltonian structures) are said to be compatible if any their linear combination 
is again a non-local Poisson vertex algebra structure (resp. a non-local Hamiltonian 
structure). Such a pair is called a bi- Hamiltonian structure. More generally, a collection 
of non-local Hamiltonian structures {H"}^^^, on V, is called compatible if any their 
(finite) linear combination is a non-local Hamiltonian structure on V. 

Recalling the Jacobi identity (4.4), we introduce the following notation. Given ra- 
tional i X ^-matrix pseudodifferential operators K,H e M.SitixiV{d), we let J{H,K) = 
J\H,K) - J^{H,K) - J''{H,K), where J^{H,K) = ( /^^.(i/, K)(A, //))^ for a = 
1, 2, 3, are the arrays with the following entries in Va,^: 

J^{H,K)ijk{X,n) = {uix{uj^Uk}H}K , 
(5.1) J'^{H,K)ijk{X,n) = {uj^{uixUk}H}K , 

J'"{H,K)ijk{X,n) = {{uixUj}Hx+^Uk}K ■ 

Consider a collection {-f^^logyt of skewadjoint rational non-local matrix pseudodifferential 
operators. By definition, H°' is a Hamiltonian structure if and only if J{H" , H°') = 0. 



22 



It is easy to see that the i?" 's form a compatible family of Hamiltonian structures if and 
only if 

(5.2) J{H'^,hP) + J{H^,H'^) = {) , \la,peA. 

Theorem 5.1. Let H, K e Mat^xf V((5) be compatible non-local Hamiltonian structures 
on the the algebra of differential functions V, which is a domain. Assume that K is an 
invertible element of the algebra Mat^x£V((?). Then the following sequence of rational 
matrix pseudodifferential operators with coefficients in V; 

i?M=ir, := (i? oK-i)""^ oi? e Mat£x£V(a) , n^l, 

form a compatible family of non-local Hamiltonian structures on V. 

Remark 5.2. It is stated in [FF81] that i^M, n ^ 0, are non-local Hamiltonian structures, 
but the prove there is given only under the additional assupmtion that H is invertible as 
well. In this case the proof becomes much easier since //['"^ is invertible, therefore one 
needs to prove that (//M)^^ is a symplectic structure. 

Following the idea in [TTll], we will reduce the proof of Theorem 5.1 to the following 
special case of it: 

Lemma 5.3. Let H, K e Mat£x£V(c5) be compatible non-local Hamiltonian structures 
on V, and assume that K is an invertible element of the algebra Mat^x£ V(5). Then the 
rational matrix pseudodifferential operator with coefficients in V 

H{d)oK-^{d)oH{d) eMat£x£V(5), 

is a non-local Hamiltonian structure on V. 

Proof. To simplify notation, in this proof we denote H by H, and we let i? = if o 

so that i?* = K-^ o H. Let iJP] = h o K'^ o h(^ = Ro H = if o i?*^, and let 

{•A'}2 = {•A'IhP] be the non-local A-bracket on V associated to ii^^ e Mat£x£V(5) 
via (4.2). We need to prove the Jacobi identity, i.e. using the notation in (5.1), that 

J(ii[2],ii[2]) =0. 

We need to compute all three terms J° = J"(iiP], iiP])jjfc(A, /x), for a = 1,2,3, of 
the Jacobi identity. First, if / e V and i e i, we have, in V((A^-'^)), 

(5.3) K;,/}2 = 2KA+a/k^i?:.(A), 

sel 

(5.4) {U,J}2 = J]{ut,+8f}H^RUl^), 

tel 

(5.5) {f\+iiUk}2 = 2 -Rfcr(A + /i + 5) {f\+^Ur}H ■ 

rel 

Both the above equations follow immediately from the Master formula (4.2) and the 
definition of if . The following identities are proved in a similar way, using that K o 
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(5.6) 
(5.7) 
(5.8) 



sel 
tel 

{fx+fiUk}H = ^Rkr{>^ + + S){fx+^Ur}K ■ 



rel 



Next, it is not hard to chek, using the left and right Leibniz rules and Lemma 3.9, that 
given an admissible non-local A-bracket {• ;\ •} on V, the following identities hold in Va,^ 
for every i,j,kel: 



(5.9) 



tel 

- ^ Rkr{>^ + fJ' + S) {uix{utyUr}K 
r,tel 

+ X Rkr{>^ + H + d) {uix{uj^Ur}H} , 



y=tJ.+d 



RtM) 



RUf^) 



rel 



(5.10) 



u 



'J /I 



sel 



^RtiW 



(5.11) 



- X i?fcr(A + /i + 5) {Uj ^{Us:^Ur} K 
r,sel 

+ X ^fcr(A + IJ. + d) {Uj^{UixUr}H} , 
rel 

{Hlf{X)x+,^Uk} = Yiii''sxUj}Hx+f,+BUk} 



x=X+d 



x=X+d 



sel 



R 

x=X+d 



:.(A)) 



- 2 {{nsxnt}Kx^,^BnkU[l^^RU^)) {\ RUf^l 

s,tel 

+ Yi^{UixUt}Hx+t,+dUk}^Rtj{f^) ■ 



tel 



Here and further we use the following notation: given an element 



N 



P(A, fi)= 2 P-,n,pA™M"(A + fif e Va^ 



m,n,p——co 



and f,geV, we let 
P{x,y) 



(5.12) 



c^X+d 



N 



2 Pm,nA^ + ^ + dr{{X + erf) {{fi + erg) e Va,, 
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In equation (5.11) we used the assumption that H and K are skewadjoint. Combining 
equations (5.3) and (5.9), equations (5.4) and (5.10), and equations (5.5) and (5.11), we 
get, respectively, 

(5.13) = {Uix{uj ^Uk}2}2 



y.{Usx{utyUk}H}H( ^ RtiWji RtjifJ-) 



s,tel 



V RkriX + fl + d) {Us^{utyUr}K}H( i?* (A) ) ( 



r,s,tel 



Yj Rkr{^ + + {Us^{Uj^Ur}H]H ( ^^^^R*i{\) 



r,sel 



(5.14) j2 _ {uj ^{u^xUk}2}2 



,tel 



y.{ujJUsxUk}H}H[ ^ RtiWji Rtjif^) 

" \ x — \ + d ' ^ 'ii—ii4-rl 



- V Rkr{\ + ^J' + ^){uty{Us^Ur}K}Hi ^ MiW] ( Rtjif^) 
r,s,tel 

+ ^kri>^ + fJ' + S){uty{uixUr}H}H( _ Rtjil^) 



r,tel 



y=H+o 



(5.15) = {{UixUj}2x+r,Uk}2 



2 -Rfer(A + H + d) {{UsxUj}Hx+^+sUr}H_ 
r,sel 

^ Rkr{>^ + H + d){{Us^Ut}Kx+^+d'<^r}H^ 



x=X+d 



r,s,tel 



x=X+d 



RU\)]Rtj{ii) 



2 Rkr{X + H + d) {{uixUt}Hx+f,+e'^r}H^Rtj{fJ-) ■ 



r,tel 



We need to prove that - - J'-^ = 0. The first term of the RHS of (5.13) combined 
with the first term of the RHS of (5.14) gives, by the Jacobi identity for H and by 
equation (5.8), 



XI [{'^sx{utyUk}H}H - {Ujy{UsxUk}H}H 



s,tel 



x=X+d 



R*M){ .RtM 



(5.16) 



y.{{UsxUt}Hnc + yUk}H}[ ^ RtiW 



sdel 



RtM 



Y,Rkr{X+^^ + S){{UsxUt}H,+yUr}K( i?*(A)Y ,^*j(/^))- 



r,s,te/ 



Similarly, the third term of the RHS of (5.13) combined with the first term of the RHS 
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of (5.15) gives, by the Jacobi identity for H and by equation (5.7), 

2 i?fer(A + l-L + d) ({Usx{uj^Ur}H}H - {{Us^Uj] H x+^.+d'^r) H 



r,sel 



(5.17) 



MiW] = T\ Rkr[>^ + H + d){Uj {UsxUr]H]H 
r,sel 

V RkriX + H + d){uty{Us.j:Ur}H}K( RtiW 



r,s.tel 



y=lJ,+d 



KM 
KM 



In the same way, the third term of the RHS of (5.14) combined with the third term of 
the RHS of (5.15) gives, by the Jacobi identity for H and by equation (5.6), 



(5.18) 



^ Rkr{^ + fl + d) y{uty{uixUr}H}H + {{Ui\Ut} H x+yUr} H 
■,tel 

(\ RtjifJ')) = - y, Rkr{>^ + fJ' + d) {uix{utyUr}H}H 

r,tel 

RtjifJ')) = - X! ^kr{>^ + fJ- + d){Usx{utyUr}H}K 



r,s,tel 



( Rti{X))( RtM 



Finally, combining the second term in the RHS of (5.13), (5.14) and (5.15), together with 
the RHS of equations (5.16), (5.17) and (5.18), we get 

- - = ^ Rkri^^ d){^- {Usx{U'tyUr^K)H 

r,s,tel 

+ {uty{UsxUr}K}H + {{UsxUt} K x+yUr} H ^ + {{UsxUt} H x+y^r} K 



+ {uty{UsxUr}H}K - {Usx{utyUr}H}K] ( RtiW 

/ \ x=X+d 

which is zero since, by assumption, H and K are compatible. 



y=H+d 



RtM 



□ 



Remark 5.4. The proof of Lemma 5.3 does not use the assumption that K \s a, Hamilto- 
nian structure. 

Proof of Theorem 5.1. We prove, by induction on n ^ 1, that the rational matrix pseu- 
dodifferential operators 

//M = K, = H,..., e MatexeV{d) , 

form a compatible family of non-local Hamiltonian structures on V. For n = 1, this 
holds by assumption. Assuming by induction that the statement holds for n ^ 1, we will 
prove that it holds for n + 1 . Namely, thanks to the observations at the beginning of the 
section, we need to prove that 

(i) J(/7["+^], //["+!]) = 0, 

(h) J(i?[™], i^["+i]) + J(i?["+i], i^M) = for every m = 0, . . . , n. 
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By the inductive assumption, H = Xir^o ^iH^^^ is a Hamiltonian structure for every 
xo, . . . ,Xn e F. Hence, by Lemma 5.3, we get the fohowing Hamiltonian structure for 
every point (xq, . . . , x„) e F"+^: 

n 2n 

where, for p = 0, . . . , 2n, 



(5.19) Qp(j;o, . . . ,x„) = ^ XiXj. 

We thus get 



ii+j=p) 



2n 



= J{HK-'H,HK-'H) = 2Qp(a;o,...,x„)J(i?W,i/W) 



p=0 

2n. 

+ 2 Qp(xo,...,x„)Q,(xo,...,x„)(J(i/W,i/M) + J(i/M,i7W)), 

{p<g) 

for every (xq, . . . , x„) e F"+-^. Note that, by the inductive assumption, J{H^\H^^) = 
for every ^ p ^ n and J(i?W, + J(i?M, = for every ^ p < q ^ n. Hence 
the above equation gives 

p=n+l 

n 2n 

(5.20) +S S Qp(^o,...,x„)Q,(xo,...,x„)(j(FM,i/M) + '/(i^f'^^f"])) 

p=0 (jr=n + l 
2n 

+ 2 Qp(2;o,...,x„)Q,(xo,...,x„)(J(/7M,/7M) + '/(^f''^,^f^^)) =0 

for every (xq, . . . , x„) e F"+^. Next, we introduce a grading in the algebra of polynomials 
in Xo, . . . , x„, letting deg(xj) = i. Then Qp{xQ, . . . , x„) is homogeneous of degree p. By 
looking at the terms of degree d = 2?i + 2 in equation (5.20), we get 

n 

Ql+l{xo, . . . , X„) J(F[-+1] , //[-+!]) + 2 Qpixo, ...,Xn) 
Q2n+2~p{xo, . . . ,X„)(j(ifM,i/P"+2-ri) + J(i/ , //M )) =0, 

while, by looking at the terms of degree d = m + n + 1 with m e {0, . . . , n} in equation 
(5.20), we get 

m 

, , /, Qp{xQ, . . . , Xn)(5m+ri+l— p(2^0i • • • 1 Xn) 

(5-22) 

(j(ifW,i?['"+"+i-rf) + J(if i/W)) =0, 
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for every {xq, . . . , Xn) £ F"^"*^. To conclude the proof, we only need to show that equations 
(5.21) and (5.22) imply respectively relations (i) and (ii) above. This is a consequence of 
the following lemma. 

Lemma 5.5. (a) For every 1, 

(5.23) Q^+i(xo, . . . , X„) ^ SpanF|(5p(xo, . . . , Xn)Q2n+2~p{xQ, . . . , x„) 

(h) For every n ^ 1 and m e {0, . . . , n}, 

(5.24) QmQn+i i Spanp] QpQ m+n+l— p 



OsSpsSm-l 



Proof. Note that, 



(5p(0, . . . , 0, x^ , . . . , Xfi) 

f if p < 2/c , 

= ^ XjXj = < x| if p = 2/c , 

«J=A: [ 2XfcXp_fc + . . . if p > 2/c . 

(i+j=p) 

We prove part (a) separately in the cases when n is even and odd. li n = 2k — 1 
is odd, letting xq = • • • = Xfc_i = we have Qn+i = x| 7^ 0, and Qp = for all p = 
2,...,n = 2k — l. This implies (5.23) for odd n. If n = 2, we have Q2 = 2xoX2 + xf, Q3 = 
2xiX2, Qi = X2, hence Q| ^ FQ2Q4- If n = 2A; with k ^ 2, letting xq = • • • = Xfc_i = we 
have <5p = for allp = 2, . . . ,n - 1, Q„ = x^, Qn+i = 2x^x^+1, Q„+2 = 2x^x^+2 + x^_^^. 
Since Q^+i = 4:x\x\_^^ is not a multiple of QnQn+2 = 2x^x^+2 + x\x\_^^, (5.23) holds for 
even Ji. 

Similarly, we prove part (b) separately in the cases when m is even and odd. If 
m = 2/c is even, letting xq = • • • = x^-i = we have QmQn+i = xf.Qn+i ¥= 0, and 
Qp = for all p = 2, . . . , m — 1. Hence (5.24) holds for even m. For m = 1 ^ n, we 
have Qo = Xg, Qi = 2xoXi, Qn+i = 2xix„ + . . . . Therefore, QoQn+2 is divisible by Xq, 
while QiQn+i = 2xoXi(2xix„ + . . . ) is not. Finally, if m = 2A; + 1 is odd, with k ^ 1, 
letting xo = • • • = Xk-i = we have Qp = for all p = 2, ... ,m — 2, Qm-i = x1, 

Qra = 2XfeXfc+i, Qn+l = 2XfcX„+i_fc + 2Xfc+iX.„_fc + . . . . HenCC, Qm-lQn+2 = xlQn+2 is 

divisible by x|, while QmQn+i = 4x|xfc+ix„+i_fe + 4xfcx|^-^x„_fc + . . . is not, proving 
(5.24) for odd m. □ 



□ 



Example 5.6. Let K = d^, H = d"^ o o ^d"^. These are compatible Hamiltonian 
structures (see [DSKWIO]). Hence, by Theorem 5.1, 



are compatible Hamiltonian structures. This was proved in [DSKWIO] by direct verifi- 
cation, and deduced from Theorem 5.1 in [TTll]. 



28 



6 Symplectic structures and Dirac structures in terms of 
non-local Hamiltonian structures 

6.1 Simplectic structure as inverse of a non-local Hamiltonian structure 

As in the previous sections, let V be an algebra of differential functions in the variables 
ui, . . . ,Ui, which is a domain, and let /C be its field of fractions. 

Recall that (see e.g. [BDSK09]) a symplectic structure on V is an ^ x ^ matrix 
differential operator S = {Sijid))- e Mat£xfV[5] which is skewadjoint and satisfies 
the following symplectic identity: 

(6.1) V f^^^A" - + (-A - M - dr^-^^) = . 

We can write the symplectic identity (6.1) in terms of the Beltrami X-bracket {■ x-} 
V X V ^ V[X], introduced in [BDSK09]. It is defined as the symmetric A-bracket such 
that (uixUj) = 6ij, and extended by the Master Formula (4.2): 



where, recalling (4.2), we let {uj-^^Uijs = Sij{X). 

Note that, if 5 e MatixeV(d) is a rational matrix pseudo differential operator with 
coefficients in V, then, by Corollary 3.11, all three terms in the LHS of equation (6.2) lie 
in Vx^fj,. Hence, equation (6.2) still makes sense (as an equation in Va,/^)- 

Definition 6.1. A non-local symplectic structure on V is a skewadjoint rational ma- 
trix pseudodifferential operator S = {Sijid)) ■ e Mat£x£V{d) with coefficients in V, 
satisftying equation (6.2) in Va,^ for all i,j, k e I. 

Theorem 6.2. Let S e Mat^x^ V(5) be a skewadjoint rational matrix pseudodifferential 
operator with coefficients in the algebra of differential functions V. Assume that S is an 
invertible element of the algebra Mat^xf V(5). Then, S is a non-local symplectic structure 
on V if and only if is a non-local Hamiltonian structure on V. 

Proof. Clearly, S is skewadjoint if and only if is skewadjoint. Hence, recalling the 
Definition 4.2 of non-local Hamiltonian structure, we only need to show that equation 
(6.1) in Va,^ is equivalent to the Jacobi identity (4.4), again in Va,^, for H = S^^. By 
equation (3.6), Remark 3.12, and the Master Formula (4.2), we have, letting Sij{d) = 




Then, the symplectic identity (6.1) becomes 



(6.2) 



(uj^{uif,Uk}s) - (uii,{uj^uk}s} + ({uj^ui}s)^^^uk) = 
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{uix{uj^Uk}H}H = {uix{S ^)kj{^J)}s-^ = 
t N 

-2 2 iS-\r{\ + f^ + d){n,,Srt,p}s-.ifi + d)PiS~%iti) 



r,t-l p=— GC> 

^^■^^ =- 2 iS-')kr{>^ + fi + d)(iX + dnS-')si{X) 



r,s,tel, neZ^ 



Exchanging i with j and A with fi, we get 

((/. + ar(5-i)t,(/.)] 



(6.4) 



\s,tel, neX-i 



du 



(n) 



Similarly, by equation (3.7) and Remark 3.12, we have, using the assumption that S is 
skewadjoint, 



{{Ui\Uj}Hx+^tUk]H = {{S ^)ji(A)A+^M/t}5-i 

t. N 



(p.K\ s,i=lp=-oc 

= 2 {s-%r{\ + ^^ + ^){-\-^x-^r 



r.s,tel, meZj 



dSts{X + d) 



(m) 



Combining equations (6.3), (6.4) and (6.5), we get that the LHS of the Jacobi identity 
(4.4) is 



{uix{Uj^Uk}H}H - {Uj^{UixUk}H}H " {{UixUj} H x+^^k} H 



(6.6) 



2 {s-%r{x + ^^ + ^)( 

r,s,teI,neZ^ \ 



du 



(n) 



du 



(n) 



where we used the notation introduced in (5.12). Clearly, the RHS of (6.6) is zero, 
provided that the symplectic identity (6.1) holds. For the opposite implication, we have, 
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by (6.6), 



2 S^fk{x + y + d)[ {Ui^{UjyUk}H}H - {Ujy{Ui^Uk}H}H 
i,j,kel \ 



-{{UixUj}Hx+yUk}H 



'~X+0 



du 



(n) 



du. 



(n) 



Hence, equation (4.4) implies equation (6.1) 



□ 



6.2 Dirac structure in terms of non-local Hamiltonian structure 

Let V be an algebra of differential functions, which is a domain, and let fC be its field of 
fractions. Given a set J and an element X e V"^ , we define the Frechet derivative of X 
as the differential operator Dx{d) : V"' given by 

neZ+ jel OUj 

Its adjoint operator is the map D%^{d) : V®'^ V®^ given by: 
(6-8) {DUd)Y)^ - 2 n-^ri^Y,). 

Here and further, for a possibly infinite set J, V®"^ denotes the space of column vectors 
in V'^ with only finitely many non-zero entries. (Though in this paper we do not consider 
any example with infinite £, we still distinguish and V®^ as a book-keeping device). 
The following identity can be checked directly and it will be useful later: 

(6.9) IX ■ Dy{d)P + ■ Dx{d)P = SP-^{X.Y), 

for all X eV-^,Ye V®-\ P e V^. 

We have the usual pairing V®^ x ^ V/dV given by {F\P) = ■ P. This pairing 
is non-degenerate (see e.g. [BDSK09, Prop. 1.3(a)]). We extend it to a non-degenerate 
symmetric bilinear form 

(6.10) <•!•>: (v®^ e vO X (v®^ ® vO ^ v/dv , 

given by <F©P|G®Q> = ^{F-Q + G-P). 

The Courant-Dorfman product is the following product on the space V®^ © V^: 

(6.11) (F © P) o (G © Q) = {DGid)P + D%{d)G - DF{d)Q + D^(a)Q) © [P, Q] , 
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where, for P,Q eV^, we let 



(6.12) [P,Q] = DQ{d)P-Dp{d)Q. 

The above formula takes a simpler form when F and G are variational derivatives (see 
[BDSK09, Rem.4.6]): 



(6.13) 




Remark 6.3. All the above notions have a natural interpretation from the point of view 
of variational calculus. Indeed, the space is naturally identified with the Lie algebra 
of evolutionary vector fields g^, and the space V®^ is naturally identified with the space 
of variational 1-forms Cl^. Then the contraction of variational 1-forms by evolutionary 
vector fields gives the inner product (6.10); the Courant-Dorfman product corresponds 
to the derived bracket [• , -J^, where [• , •] is the Lie superalgebra bracket on the space of 
endomorphisms of the space of all de Rham forms over V, and d = ad((5), where 5 is the 
de Rham differential, [BDSK09, Prop.4.2]. 

Definition 6.4 ([Dor93, BDSK09]). A Dirac structure is a subspace £ c V®*^ © 
which is maximal isotropic with respect to the inner product (6.10), and which is closed 
under the Courant-Dorfman product (6.11). 

Given two ixi matrix differential operators A,Be Matexi consider the following 
subspace of V®^ © V^: 

(6.14) Ca,b = {B{d)X®A{d)X\X eV®^}. 

Proposition 6.5. The subspace Ca,b V®^ © is isotropic with respect to the inner 
product (6.10) if and only if 

(6.15) A*oB + B*oA = 0. 

If, moreover, deti? 7^ (i.e. B is invertible in the algebra Mat£xeK^{{d^^))), then (6.15) 
holds if and only if Ao B^^ e Mati^e is skewadjiont, while if det A ^ 0, then 

(6.15) holds if and only if B o A^^ e Mat^xe is skewadjiont. 

Proof For X,Y e V®^ we have 

(B{d)xeA{d)x\B{d)Y®A{d)Y) = • {A*{d)Bid) + B*{d)A{d))x . 

Hence, due to non-degeeracy of the pairing {F\P) = S^F ■ P, the space Ca,b is isotropic 
if and only if (6.15) holds. The remaining statements are straightforward. □ 

Example 6.6. Letting A e Mat^ V and B = l^d, condition (6.15) holds if and only if A is 
a symmetric matrix with entries in C V (the subring of constant functions) . In this case 
AB~^ is a skewadjont matrix pseudodifferential operator and Ca,b = {AX © dX | X e 
V®^} is an isotropic subspace of V®^ © V^. It is not hard to show directly that Ca,b 
is maximal isotropic if and only if the matrix A is non-degenerate. When V = /C is a 
differential field, this is a corollary of the following general result: 
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Proposition 6.7 ([CDSK12b]). Let K he a differential field, and let H = AB^^ be a 
minimal fractional decomposition of the skewadjoint rational matrix pseudodifferential 
operator H e Mat^x£^(f^)- Then the subspace Ca,b ^ /C®^ © /C^ is maximal isotropic 
with respect to the inner product (6.10). 

Proposition 6.8. Suppose that A,B e Mat^x^V[5] satisfy equation (6.15). Then the 
following conditions are equivalent: 

(i) <X o y, Z> = for all X,Y,Ze Ca,b- 

(a) for every F,G eV^ one has: 

A*{d)Ds^S)Gid)A{d)F + A%d)Dl^,^p{d)B{d)G 
(6.16) -A-{d)DB(8)F{.S)A{d)G + A-{d)Dl^g^^{d)A{d)G 

+B*{d)DA^S)G{d)A{d)F - B*{d)DAi^8)F{d)A{d)G = 0. 

(Hi) for every i,j,k e I, one has in the space V[A,/i].' 

S f^L(A + /X + 5) f ^^(A + erAuix) - ^^(/^ + drA.if^)) 

neZ+ 

(6-17) +BliX + + a) f ^^(A + dTAuiX) - ^M^(^ + 5)M,,(/.)) 

\ dul oui J 

+AUx + ^ + ^){-X-^,-^r(^^B,,{^,) + ^-^A,,{^,)^ = o. 

V ous dus ' J 

Proof. Letting X = B{d)F@A{d)F, Y = B{d)G@A{d)G, Z = B{d)E®A{d)E, condition 
(i) reads 

l{A{d)E) ■ {DBid)G{d)A{d)F + D\^^p{d)B{d)G - DB(^8)F[d)A{d)G 
+D*(5)^(5)^(5)G) + (i?(5)ii;)-(l)^(5)G(5)^(5)F-Z?^(,)^(a)A(a)G)=0. 

Since the above equation holds for every E e V®^ it reduces, integrating by parts, to 
equation (6.16). For this we use the non-degeneracy of the pairing (6.10). This proves 
that conditions (i) and (ii) are equivalent. 

We next prove that conditions (ii) and (iii) are equivalent, provided that (6.15) 
holds. For a = 1,...,6, let (6.16)q, be the /c-entry of the a-th term of the LHS of 
(6.16): for example (6.16)i = {A* {d)DB[S)G{S)A{d)F) ^. We have, by the definition of 
the Frechet derivative and some algebraic manipulations (similar to those used in the 
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proof of [BDSK09, Prop.1.16]), 



(6.16)i =22 (AtM{^^G,)d-Au{d)F, 

i,j,s,tel neZ+ ^ ^U^ 
r)C ■ 

+AUd)Bs,{d)^^d-Au{d)F, 

(6.16)2 =22 {AlA^){-^T{^-^F^|B,,{^)G, 

dFi 



i,j,s,tel neX 



du 



+AtM{-dr-j^At,{d)Bt,{d)G, 

OUs 

(6-16)3 = -22 (AtsiS){^^F)d-A,id)G, 

i,j,s,tel neZ+ ^ du^ 
dF- 

+AUd)Bs^{d)^d^At,{d)G, 

du\ 

(6.16)4 =22 iAlsm-dr{^^F)A,,{d)G, 

i,j,s,tel neZ+ ^ ^^^s 
dF 

+AUd){-dy^^BmAt,{d)G, 
dus 

(6.16)5 =22 iBUd){^^G,)d-Au{d)F, 

i,j,s,tel neZ+ ^ du^ 



dG ■ 

+BUd)A,,{d)^d-Au{d)F, 



(6.16)6 = -22 (BW){^f^F,)d-AMG. 

i,j,s,tel neZ+ ^ C^"^* 
dF- 

+Bl{d)AAd)^d-At,{d)G, 

du\ 

Combining the second terms in the RHS of (6.16)i and (6.16)5 we get zero, thanks to 
equation (6.15). Similarly, we get zero if we combine the second terms in the RHS of 
(6.16)2 and (6.16)4, and if we combine the second terms in the RHS of (6.16)3 and (6.16)6- 
Equation (6.17) thus follows from equation (6.16) once we replace d acting on Fi by A, 
and d acting on Gj by /i. □ 

Remark 6.9. It follows by Definition 6.4 and Proposition 6.8 that a Dirac structure is 
a maximal isotropic subspace C of V®^ © satisfying one of the equivalent conditions 
(i)-(iii) above. 

Proposition 6.10. Suppose that A, B e Mat£x£V[5] satisfy equation (6.15) and assume 
that B has non-zero Dieudonne determinant. Suppose, moreover, that the (skew adjoint) 
rational matrix pseudodifferential operator H = Ao B~^ has coefficients in V, i.e. He 
Mat£x£ V(5). Consider the corresponding non-local X-bracket {■ \ given by the Master 
Formula (4.2). Then the Jacobi identity (4.4) on {• \ ■}h is equivalent to equation (6.17) 
on the entries of matrices A and B. 
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Proof. Letting Ast{d) = J^^^o "■st;md'^ and Bst{d) = Xi^lo &ji;n^'"- By formula (4.3) and 
the left Leibniz rule (3.4) we have, 

{uix{Uj^Uk}H}H = 2{^iiA^fcr(M + d) 3:^^^ (fl)} H 
rel 

M 

= 2 S {^iAafcr-;m}H(/^ + d)"'B;^{fl) 
rel m^O 

+ 2 Akr{X + + d){u^xB;l{fi)}H ■ 

rel 

By equation (3.8) we have 

M 

= - S 2 (B-%s{X + ^ + ^){u,xbst■,r,^}Hi^^ + dy^iB-Xijl) . 
s,tel m=0 

Combining the above two equations we then get, using the Master Formula (4.2), 



H 



M 



2 2 T.^-^{i^^ + srB;li^.))[ix + ^rA,,ix + ^)B-^^ 

r,s,teIm=OneZ+ SUg 



M 



S H 2 Mr{X + H + d){B~\s{X + ^i + d) 



p,q,r,s,tel m=0 neZ+ 



db 



^((^ + dr{B-%{f,)) ((A + drA,,{X + d)B;^{X)) 
2 2 C^'iit ^^ ^n (^)) ((^ + SrAAX + d)B,HX)) 

r,s,teIneZ+ (J^s 

- S Mr{X + 11 + d){B-\s{X + ^x + d) 

p,q,r,s,tel neX-f- 



dBst{fi + d) 



iB-%i^,)) ((A + drA„{X + d)B;^iX)) 



Next, we apply B1^,j^{X + + c5) to both sides of the above equation (on the left), replace 
A by A + 5 acting on Bai^X), replace fj, hy fi + d acting on Bjji{p), and sum over i,j,kel. 
As a result we get, using the assumption (6.15) (see notation (5.12)), 

y! Bl,i.{X + fi + d){ui^{uj Uk}H}H( ^ Bii>{X))( Bjjfii^i) 



i,j,kel 



(6.18) = 2 S {Bt,k{X + H + d)^-^^{X + drAu'{X) 

+Al,^{X + ^ + + drAu>{X) 
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Exchanging i' with j' and A with ^ in (6.18), we get 

2 Bl>k{^ + + S){Ujy{UixUk}H}H 
i,j,kel 

= 2 ^k'ki^ + + S){Uix{UjyUk}H}H 
i,j,kel 



Bii'{\) 



Bij'ifJ.) 



Bji'{X) 



(6.19) 



+A^,,(A + /. + a)^^^(M + r?)M,y(/i) 

We are left to study the third term in the Jacobi identity. By the right Leibniz rule (3.4) 
we get, 

{{UixUj}Hx+^Uk}H = Yj{Ajri\ + d)B~^^ {X) x+f,Uk} H 



rel 



M 



S S {^^r;m,^,+8^k}HjX + d)"^B;^{X) 

rel m=0 
M 

2 2 {B;^{X)x+^,+eUk}H^{-l^ - drajr;r 



rel m=0 

By equation (3.9) we have 

{Bri^ {X)x+^+0Uk} H ^ 
I M 

s,t=l m=0 

Combining the above two equations and using the Master Formula (4.2) we then get 

{{uixUj}Hx+^Uk}H 



2 Aktix + /X + d)Br,\x + ^ + d){-x - ^ - dr^^^^^^ + 



r,s,tel 



(5n, 



(n) 



- 2 ^fc,(A + /i + a)i?-pi(A + // + 5)(-A-//-5)'^ 



r,p,q,s,tel 
neZ+ 



5ti 



^(/X + 5)^*-(/i) 



Hence, if we apply, as before, B'^,j^{X + fi + d) on the left, replace A by A + 5 acting on 
Biii{X), replace ^hy ^ + d acting on Bjji{fi), and sum over k e I, we get, using (6.15), 



2 Bli.{X + fi + d){{ui^Uj}H^+yUk}H 

i,j,kel 



x=X+d 



Bii> (A) 



(6.20) 



Bjj'il^) 



- 2 2 Al,,{x + ^, + ^){-x-^,-^r 

j,kel neZ-i- 
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Combining (6.18), (6.19) and (6.20), we get that the expression 



^k'ki^ + fJ- + S)({ui^{ujyUk}H}H - {ujy{Ui^Uk}H}H 



i,j,kel 



is the same as the LHS of (6.17). The claim follows from the assumption that the matrix 
B e Mat^x£V[(?] has non-zero Dieudonne determinant. □ 

Theorem 6.11. Let V be an algebra of differential functions in the variables ui, . . . , ui, 
which is a domain, and let JC be its field of fractions. Let H = AB~^ , with A, B e 
M.at£xeV[d], detB ¥= 0, be a minimal fractional decomposition (cf. Definition 2.11 and 
Remark 2.14) of the rational matrix pseudodifferential operator H e Mat^x^ V(5). Then 
the subspace 

(6.21) CaA^) = {B{d)X ® A{d)A \ X e /C®^} c /C®^' ® /C^ , 

is a Dirac structure if and only if H is a non-local Hamiltonian structure on V. 

Proof. It immediately follows from Remark 2.14 and Propositions 6.5, 6.7, 6.8 and 6.10. 

□ 

Remark 6.12. We may define a "generalized" Dirac structure subspace/: of V®^©V^ 
such that C c (i.e. C is isotropic), and Co C a (i.e. condition (i) in Proposition 
6.8 holds), where is the orthogonal complement to L with respect to the inner product 
(6.10). Note that a Dirac structure is a special case of this when £ is maximal isotropic. 
\i A^B e Mat^x^ deti? i= 0, then £a,b is a generalized Dirac structure if and only 
\i H. = AB~^ is a non-local Hamiltonian structure on V (not necessarily in its minimal 
fractional decomposition). Note also that any subspace of a generalized Dirac structure 
is a generalized Dirac structure. 



6.3 Compatible pairs of Dirac structures 

The notion of compatibility of Dirac structures was introduced by Gelfand and Dorfman 
[GD80], [Dor93] (see also [BDSK09]). In this paper we introduce a weaker, but more 
natural, notion of compatibility, which still can be used to implement successfully the 
Lenard-Magri scheme of integrability, and which is more closely related to the notion of 
compatibility of the corresponding non-local Hamiltonian structures. 

Given two Dirac structures L and CJ cz V®^ ® V^, we define the relations 

F®PeC,F@P' e C for some F e V®^} c © V^ 
F®PeC,F'@PeC' for some P e V^} c V®^eV®^. 

Definition 6.13. Two Dirac structures C c V®'^ ® V'^ are said to be compatible if 
for all P,P',Q,Q' eV^, F,F', F" e V®^ such that 

P®P',Q®Q' e Mc,c' and F © F', F' @ F" e Mc,c' , 



(6 22) ^c,c' = {P®P' 



37 



we have 

(6.23) iF\[P,Q]) - iF'\[P,Q']) - {F'\[P',Q]) + {F"\[P',Q']) = 0, 
where, as before, {F\P) = • P, and, for P,Qe V\ [P, Q\ is given by (6.12). 

Remark 6.14. The original notion of compatibility, introduced by Dorfman in [Dor93], is 
similar, except that Mc,c' replaced by the "dual" relation 

^C,C' = {F®F' e V®^ © V®^ I • P = \F' ■ P' for all P ® P' e Mc,c'} ■ 

Since C and C are isotropic, we have, for F © F' e Mc,c'-> and for Q @ Q' e Mc,c'i 
\F -Q = -\G -P = \F' ■ Q', where P e V^' and G e V®^ are such that F © F, G ® Q e 
£, F' © F, G © Q' e C. Hence, Mc,c' ^ ^c,c'- 

Even with the weaker notion of compatibility, the following important theorem still 
holds (cf. [BDSK09, Thm.4.13]). 

Theorem 6.15. Let (£, C) he a pair of compatible Dirac structures. Let Fq, Fi, F2 e V®^ 
be such that: 

(^) D^fS^) = DfM, /orn = 0,1; 
(ii) Fo © Fi , Fi © F2 e Mc,c' ■ 
Then, for all P@P',Q@Q' e Mc,c', we have 

(6.24) IQ'-{DfM-D*fMP' = 0- 
Proof. By the assumption (6.23), we have 

= (Fo|[F,Q]) - (Fi|[F,Q']) - (Fi|[F',Q]) + (F2|[F',Q']) 
= J (Fo • DQ{d)P - Fo • Dp{d)Q - Fi • DQ,{d)P + Fi • Dp{d)Q' 

-Fi • DQ{d)P' + Fi • Dp>id)Q + F2 • DQ>{d)P' - F2 • Dp>{d)Q'] 
= \P- l^a^olQ) - {F,\Q')) - Jq • ^{{Fo\P) - {Fi\P')) 
-]P'- ^((^ilQ) - (^2|Q')) +\q'- ^mP) - (^2 IF')) 
- ^ Q ■ DfMP + ^ P ■ Dfo {S)Q + ^Q'- Dp, {d)P -\^P-Dp, {d)Q' 
+ Jq • DF,{d)P' - J^' • DF,{d)Q - Jq' • ^F2(5)F' + J^' • ^F2(5)Q' . 

In the second identity we used the definition (6.12) of the Lie bracket on V^, and in 
the last identity we used equation (6.9). Since, by assumption, Fo © Fi e Nc,c' aiid 
Q@Q' e Mc,c', we have (by Remark 6.14) that (FqIQ) = (Fi|Q'). Hence the first term 
in the RHS above is zero, and, by the same argument, the first four terms are zero. The 
following six terms are also zero since, by assumption, DFg(d) and DF-^{d) are selfadjoint. 
In conclusion, equation (6.24) holds. □ 
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6.4 Compatible non-local Hamiltonian structures and corresponding 
compatible pairs of Dirac structures 

In Theorem 6.11 we proved that to a non-local Hamiltonian structure H e Mat^x^ V((5) 
in its minimal fractional decomposition H = AB~^, with A,Be Mat^x^V[5], there 
corresponds a Dirac structure CA,B{Ki) on the field of frations /C. In this section we prove 
that to a compatible pair of non-local Hamiltonian structures H = AB~^, K = CD~^, 
in their minimal fractional decompositions, there corresponds a compatible pair of Dirac 
structures Ca,b{IC), J^c,d{K,) on /C. This is stated in the following: 

Theorem 6.16. Let V be an algebra of differential functions in ui, . . . ,ui, which is a 
domain, and let IC be its field of fractions. LetH, K e Mat^x^ V(5) be compatible non-local 
Hamiltonian structures on V. Let H = AB~^^ K = CD^^ be their minimal fractional 
decompositions (cf. Definition 2.11). Then Ca,b{K^) and Cc,d{K,) are compatible Dirac 
structures on IC. 

By Theorem 4.1, the Hamiltonian structures H and K on V are compatible if and 
only if we have the following "mixed" Jacobi identity on generators {i,j, k e L): 



In order to relate the above condition to the compatibility of the corresponding Dirac 
structures Ca,b and Cc,d-, we need to compute explicitly each term of the above equation. 
This is done in the following: 

Lemma 6.17. Suppose that the pairs {A, B) and (C, D), with A, B,C,D e Mat^x^ "^'[^1; 
satisfy equation (6.15): 



Assume that B and D have non-zero Dieudonne determinant, and that the (skewadjoint) 
rational matrix pseudodifferential operators H = AB^^ and K = CD~^ have coefficients 
in V, i.e. H,K e Mat^xf V((?). Consider the corresponding non-local X-brackets {■ \ ■}h 
and {-x-}!^ given by the Master Formula (4.2). Then, in terms of notation (5.12), we 
have the following identities for every i' ,j' , k' e L: 



(6.25) 



{Uix{Uj^Uk}H}K - {Uj^{UixUk}H}K " {{UixUj}H x+^'^^kjx 

+ {'>J'ix{Uj^Uk}K}H - {Uj^{UixUk}K}H " {{UixUj}K x+^Uk}H = 



(6.26) 



A* 



oB + B*oA = 0, C*oD + D*oC = 0. 



2 Bk'ki>' + + ^){uix{ujyUk}H}K{ 



x=X+d 



aA^'(A))( 



Bjfifi)) 



i,j,kel 



= 22 Bl,,iX + f, + d) 



(6.27) 



i,kel neZ+ 
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(6.28) 



(6.29) 



(6.30) 



(m)) 

i,j,kel 

= S S Dt,,iX + ^^ + d)^-^^{X + orA,AX) 

i,keIneZ+ ^U- 

i,kel neZ+ ^U- 

Yi ^k'ki^ + + S){ujy{ui^Uk]H]K{\^^y^^gB,,,{X)) (|y=^+5^ij'(^)) 
i,j,kel 

= 22 Bl,,{x + ^, + ^)^-^^[^, + ^TC,,,{^,) 

j,keIneZ+ ^Uj 

+ 2 J] At,,{x + ^^ + ^)^^^i^^ + ^rc,yi^^), 

2 + ^ + ^){%y{'"ix^ifc}i^}//(|^^;^+5Ai'(A)) (|j^^^^5-Bjy(^)) 

+ 2 2 ^.%(A + /^ + 5)^^^(/. + 5)M,,'(/x), 

j,keIneZ+ du^ 



i,j,kel 



ij,kel 

(6.31) = - S S + + 5)(-A - /. - a)"^^^i?,y(^) 

j,keIneZ+ VU^ 



j,keIneZ+ SUf, 
X! ^fc'fc(A + M + 5){Kx^^ik^+y^^fc}H(|^^;,+5Ai'(A))(|^^^^^,I)jy(/i^ 



(6.32) = - S S ^^''^(^ + + a)(-A - /. - ^y^^^^D,,,{^,) 

- S S ^^'^(^ + + 5)(-A - ^ - a)'^^^^c,H/^) . 

j,keInsZ+ SU/^ 

Proof. For equation (6.27), we can use the Leibniz rule and equation (3.8) to get 

{Uix{UjyUk}H}K = Y^{ui^Akr{y + d){B^'^)rj{y)} K 
rel 

= S S {^ixakr;m}K{y + dT{B-^)rj{y) 
rel meZ-f- 

- Yi ^kr{x + y + d){B-\.p{x + y + d){uixbpq;m}K(.y + d)'^{B-^)gj{y) . 



r,p,qel 
meZj 
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We can then use the Master Formula (4.2) to get 



{UiAujyUklHlK = 2j Zj ( Tl^) 



r,p,q,sel 



r,sel neZ+ OUg 




If we now replace x with \ + d acting on Dai (A) and yhy ^ + d acting on Bjji [i^l) , and we 
apply -B^/fc(A + ^ + d), acting from the left, to both sides of the above equation, we get, 
after using the assuption (6.26), that equation (6.27) holds. Equation (6.28) is obtained 
from (6.27) by exchanging the roles of H and K. Equation (6.29) is obtained from (6.27) 
by exchanging A with jj. and i and i' with j and j' respectively, and equation (6.30) is 
obtained from (6.29) by exchaing the roles of H and K. Finally, equations (6.31) and 
(6.32) can be derived with a similar computation, which involves the right Leibniz rule 
(instead of the left) and equation (3.9) (instead of (3.8)). □ 



Let us next describe the relations (6.22) associated to Dirac structures Ca,b{K-) and 
C.C,d{K,) defined in (6.21). We have 

Mc^ smxcoiic) = {Md)X e C{d)X' \X,X'e K®' , B{d)X = D{d)X'} , 

Hence, by Definition 6.13, the Dirac structures Ca,b and Cc^d are compatible if and only 
if, for every X, X\ Y, Y' , Z, Z' , W, W e V®^ such that 

B{d)X = D{d)X' , B{d)Y = D{d)Y' , B{d)W = D{d)Z' , 

(6.34) 

A{d)z = c{d)z' , A{d)w = cid)w' , 



Lemma 6.18. Suppose that H = AB^ and K = CD are non-local Hamiltonian 
structures, and that conditions (6.34) hold. Then equation (6.35) is equivalent to the 
following equation: 



- l{A{d)Y) . DB^s)xid)Cid)Z' + mS)Y) ■ D%^,^^id)Cid)Z' 
+ l{B{d)Y) ■ Dc^sjz'{d)A{d)X - mS)Y) ■ DAie)x{d)C{d)Z' 
+ S{C{d)Y') ■ Dnis)Z'{d)A{d)X + \{C{d)Y') • D\^,^^{d)D{d)Z' 
+ mm ■ DBis)w{d)C{d)X' + l{A{d)Y) . Dl.^g^^,{d)B{d)W 

- KC(a)y') • D^^s^x,{d)A{d)W + mS)Y') • D},^,^^,{d)A{d)W 

+ \{D{d)Y') . DAi8)wid)C{d)X' - S{D{d)Y') • Dc(s)x'{d)A{d)W = . 



we have the following identity: 



(6.35) 



{B{d)Z\[A{d)X,A{d)Y]) - {D{d)Z'\[A{d)X,C{d)Y']) 

- {B{d)w\[Cid)x',A{d)Y]) + {D{d)w'\[cid)x' ,c{d)Y']) = o. 
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Proof. By (6.12) and (6.9), we have 

{Bid)Z\[A{d)X,A{d)Y]) = \{B{d)Z) ■ DAie)Y{d)A{d)X 
(6.37) - ^{B{d)Z) ■ DAie)x{d)A{d)Y = S(A(a)X) • ^{B{d)Z\A{d)Y) 

- liA{d)Y) ■ DB^s)z{d)A{d)X - l{A{d)Y) ■ D*^^,^^{d)B{d)Z . 

Similarly, we have 

{D{d)Z'\[A{d)X,C{d)Y']) = • j^{D{d)Z'\C{d)Y') 

- KC(a)y') • D^^p),,{d)A{d)X - ^{Cid)Y') ■ D\^g^^{d)D{d)Z' , 



{B{d)W\{C{d)X\A{d)Y]) = \{C{d)X') . j-{B{d)W\A{d)Y) 
- mS)y) ■ DBi^s)w{S)C{d)X' - l{A{d)Yy D},^g)^,{d)B{d)W , 

and 

{D{d)W'\[C{d)X',C{d)Y']) = \{C{d)X') . -^{D{d)W'\C{d)Y') 
- ■ DniS)W'{d)C{d)X' - KC(a)yO . Dl^s)AS)D{d)W' . 

By the skewadnointness of H and K, which translates to (6.26), and by conditions (6.34), 
we have 

{B{d)Z\A{d)Y) = -{A{d)Z\B{d)Y) 

= -{C{d)Z'\D{d)Y') = {D{d)Z'\C{d)Y'), 

hence the first terms in the RHS of (6.37) and (6.38) cancel. Similarly for the first terms 
in the RHS of (6.39) and (6.40). Therefore, combining equations (6.37)-(6.40), we get 
that equation (6.35) is equivalent to 

- l{A{d)Y) . DB^8)z{d)A{d)X - l{A{d)Y) ■ Dl^s)xm{d)Z 

+ \{C{d)Y') . Dj,^ff)Z'{d)A{d)X + l{C{d)Y') ■ D*^^g^^{d)D{d)Z' 

(6.41) 

+ l{A{d)Y) ■ Ds(8)w{d)C{d)X' + l{A{d)Y) ■ D*^,^^,id)Bid)W 

- mS)y') ■ Dj,^s)W'{d)C{d)X' - S{C{d)Y') ■ D}.^,^^,{d)D{d)W' = 0. 

Next, since by assumption H = AB^^ is a non-local Hamiltonian structure, it follows by 
Propositions 6.10 and 6.8 that equation (6.16) holds. In particular, 

- s{A{d)Y) ■ Dj,^s)zid)Aid)x - ms)Y) ■ D\^^s)xmm 

(6.42) = -l{A{d)Y) . Dsis)x{d)A{d)Z + mS)Y) ■ D%^,^^{d)A{d)Z 

+ l{B{d)Y) . D^(^8^z{d)A{d)X - ' I?A(5)x(5)A(a)Z . 
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Similarly, using the assumption that K = CD is a non-local Hamiltonian structure, 
we get 

- \{C{d)Y') ■ Dj,^s)W'{d)C{d)X' - md)y') ■ Dl^^s)x>m{d)W' 
(6.43) = -l{C{d)Y') . DD(^s)X'{d)C{d)W' + l{C{d)Y') ■ Dl^s)X'{S)C{d)W' 

+ l{D{d)Y') ■ Dc(^s)w'{d)Cid)X' - \{D{d)Y') ■ I)c(5)X'(5)C(5)VF' . 

Combining equations (6.41), (6.42) and (6.43), we get (6.36). □ 

Proof of Theorem 6.16. By Lemma 6.18, we only need to prove that, if condition (6.25) 
holds, then equation (6.36) holds for every X,X' ,Y,Y' ,W,Z' satisfying the first three 
identities in (6.34). It follows by some straightforward computation that we can rewrite 
each term in the LHS of (6.36) as follows 

- S{A{d)Y) ■ Dsi^s)x{d)C{d)Z' = -S{A{d)Y) • B{d)Dx{d)C{d)Z' 

i',j',k'el 
j,k,el, n.eZ+ 



S{Aid)Y) ■ Dl^,^^,{d)C{d)Z' = S{A{d)Y) ■ D\,{d)D*{d)C{d)Z' 



'i',j',k'el ^"-k 
j,k,eI,neZ+ 



S{D{d)Y') ■ Dcis)z' mm = ^{Dm') ■ cm z' mm 

(6-^6) +f 2 n'^k(A + ^^ + 5)^^^(A + 5)M.,(A)(|,^,X,) (1,^,4) 



i',j',k'el 
i,k,el jneX-f. 



- sism) ■ DAi8)xmm' = -SimY) ■ Amxmm' 

(6-47) -f 2 n,i?,\(A + ^ + a)^^(/. + 5)-C,,,(/^)(L.,^.)(L.,4) 



i',j',k'el ^"-j 
j,k,eI,neZ^ 



\{cm') ■ DDi8)z'mm = licm') ■ omz' mm 

+\ 2 n''^^.%(A+^+5)^^^(A + 5)"yl.,(A)(|,_,X,)(|^_,4) 



i',j',k'el 
i,k,el ,neZ^ 
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(6.49) + r 2 n,c:4X+^^+^){-x-^.-^r^-^^B,,i^^){\^^^^^ 



i',j',k'el ^'^k 
j,k,el ,neZ+ 



K^(5)y) • DB^s)wid)C{d)X' = S{A{d)Y) ■ B{d)Dw{d)C{d)X' 

i',j',k'el 0U\ 
i,k,el ,nel,^ 



mm) ■ D*c^s)X'iS)D{d)Z' = S{A{d)Y) ■ D\,{d)C%d)D{d)Z' 
•'i'J^k'el ^% 



- KC(5)y') • Dn^s)x'm{d)W = -md)y') ■ D{d)Dx'id)Aid)W 

i',j',k'el ^""j 
j,fc,G/,nGZ+ 



''i',j',k'el "^k 
j,k,el ,neZ+ 



mm) ■ DAi3)wid)C{d)X' = mm) ■ A{d)Dw{d)C{d)X' 

(6-54) +f 2 YyBl,,{X + ,. + d)^-^^{X + drCAmx-8^M,-^^^^^ 
i',j',k'si du\ 

i,k,el ,neZ^ 



- mim') ■ Dc^s)j,,{d)A{d)w = -mim') ■ c{d)Dx'{d)A{d)w 

(6-55) -f 2 yk^Dt,,ix+^.+^)^-^^i^. + ^rA,,i^^^^^ 

j,k,el ,neZ^ 

It follows from the skewadjointness conditions (6.26) that the first term in the RHS of 

(6.44) cancels with the first term in the RHS of (6.47), the first term in the RHS of 

(6.45) cancels with the first term in the RHS of (6.51), the first term in the RHS of 



44 



(6.46) cancels with the first term in the RHS of (6.48), the first term in the RHS of 

(6.49) cancels with the first term in the RHS of (6.53), the first term in the RHS of 

(6.50) cancels with the first term in the RHS of (6.54), and the first term in the RHS 
of (6.52) cancels with the first term in the RHS of (6.55). Furthermore, combining the 
second terms of the RHS's of (6.50) and (6.54), we get, thanks to (6.27), 

f 2 (i?fcfe.(5)n,){na{^x,^nfc}^,}K(L^,AiK5)X.'0(L.5^.y(^)W^/) ■ 

t,j,kel 

Combining the second terms of the RHS's of (6.46) and (6.48), we get, thanks to (6.28), 
[ 2 (^fcfc'(5)y;){niA{n,-^nfc}i^}H(L=5i?n'(5)^^')(L=a^jy('^)4) • 

i',j',k'el 
i,j,k,el 

Combining the second terms of the RHS's of (6.44) and (6.47), we get, thanks to (6.29), 

Ti {Bkk'{d)Yk>){uj^{uixUk}H}K{\^^gBii,{d)Xi,){\^^gDjj,{d)Zj,) . 

•'i',f,k'el 

i,j,k,el 

Combining the second terms of the RHS's of (6.52) and (6.55), we get, thanks to (6.30), 

- [ 2 (I)fcfc.(a)n!,){^.7.{«a^fc}ir}i/(L.5A^'(5)X,',)(|^^^i?,,,(a)l^,v). 

''i',j',k'el 
i,j,k,el 

Combining the second terms of the RHS's of (6.49) and (6.53), we get, thanks to (6.31), 

- [ 2 (l)fefc'(a)yfc',){{nan,}m+M7Xfc}x(L=5i?n'(^5)^.')(L=a%'(^)^/)- 

i',f,k'el 
i,j,k,el 

Finally, combining the second terms of the RHS's of (6.45) and (6.51), we get, thanks to 
(6.32), 

- [ 2 (i?fcfc<5)n,){{7x,;,^.,};,;^^^n,}^^(|^^^A^'(^)^.')(L=<5^.y(^)4)- 

i,j,k,el 

Putting together all the above results, we conclude that the LHS of (6.36) is equal to 

- {Uj^{uixUk}H}K - {uj^{uiXUk}K}H " {{uixUj} H x+^^k} K 

-{{u.xUj}Kx+,^k}H){\x^sBid)X)^{l^^Bid)W)^, 
which is zero by (6.25). □ 
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In view of Theorem 6.16, we can translate Theorem 6.15 in terms of compatible 
non-local Hamiltonian structmes. 

Theorem 6.19. Let V be an algebra of differential functions in ui, . . . ,U£, which is 
a domain. Let H,K e Ma.tixeV{d) be compatible non-local Hamiltonian structures 
on V. Let H = AB^^, and K = CD^^ be their minimal fractional decompositions 
(cf. Definition 2.11), with A,B,C,D e Mat^x^V[a], detB ^ 0, detZ) 0. Let 
Fq = B{d)Z, Fi = D{d)Z' = B{d)W, F2 = D{d)W', with Z,Z',W,W' e V®^ be such 
that 

D{d)Z' = B{d)W , A{d)Z = C{d)Z' , A{d)W = C{d)W' , 

and 

D%{d) = DF,{d) , D%^{d) = DF,{d). 

Then: 

(a) For all X,Y e V®^ such that D{d)X, D{d)Y e lm{B), we have 
(6.56) lY ■ C* (d) [Dp, (d) -D*p^{d))C{d)X = 0. 

(b) If we also assume that detK 7^ 0, then Dp^{d) = Dp^i^S). 

(c) If, moreover, we assume that V is a normal algebra of differential functions, then F2 
is exact: F2 = ^ for some 5/2 e V/dV. 

Proof. By Theorem 6.16, Ca,b{K^) and Cc,d{K-) are compatible Dirac structures over /C, 
the field of fractions of V. Recalling the expressions (6.33) of Mc,C' and Afc^c for these 
Dirac structures, we get by Theorem 6.15 that equation (6.56) holds over /C, hence over 
V, proving (a). Let us prove part (b). It is proved in [CDSK12b] that any two matrices 
B[d) and D{d) with non zero determinant have a right common multiple B{d)Di{d) = 
D{d)Bi{d), where Bi{d), Di{d) e Mat^x£^[5] have non-zero determinant. By clearing 
denominators, we can assume that Bi{d) and Di{d) have coefficients in V. Hence, if 
X,Y e Im(i?i), we have D[d)X,D{d)Y e \m{B). Therefore, by part (a) we have 

Jg • Bl{d)C''{d){DFM - Dl^{d))C{d)Bi{d)F = 0, 

for all F,G e V®^. Since, by assumption, G and Bi have non-zero determinants, it 
follows that Dp^(d) = DF^{d), as we wanted. Finally, part (c) follows by the fact that, 
under the assumption that V is normal, the variational complex is exact (see [BDSK09, 
Thm.3.2]). □ 

7 Hamiltonian equations associated to a non-local Hamil- 
tonian structure 

7.1 Local functionals, Hamiltonian equations and integrability 

Let V be an algebra of differential functions in the variables Ui, i e I, assume that it is a 
domain, and let /C be it field of fractions. Let H e Matixe ^{(^) be a non-local Hamiltonian 
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structure on V. Recall that, since H has rational entries, it admits a minimal fractional 
decomposition H = AB^^, where A,B e Matixt V[d] and det B (cf. Definition 2.11). 
Throughout this section we fix such a minimal fractional decomposition for H. Recall 
that Ca,b is a Dirac structure on V by Theorem 6.11. 

Definition 7.1. A Hamiltonian functional (for H = AB~^) is an element e V/dV 
such that 1^ = B{d)F for some F e V®^. Then, P = A{d)F e is called a Hamiltonian 
vector field associated to ^h. We denote by J-{H) cz V/dV the subspace of all Hamiltonian 
functionals, and by ^{{H) a the subspace of all Hamiltonian vector fields: 

T{H) = (^^y\lmB) ^V/dV , = ^(^-^(im^)) c 

We say that J/ e and P e T~i{H) are associated (with respect to the Hamiltonian 

structure H = AP-^) \i there exists F eV^ such that |{ = S(5)F, P = yl(5)F. In this 

case we write \f P. 

Remark 7.2. Note that, the spaces F{H) cz V/dV and 'H{H) a V^, as well as the 

relation \h P, for J/i e /"(i?) and Pen{H), may depend not only on the non-local 
Hamiltonian structure H, but also on its fractional decomposition H = AB~^. However, 
if A and B are multiplied on the right by a matrix differential operator D which is 
invertible in the algebra Mat^x^ V[(?], the spaces J~{H) and 'H{H), as well as the relation 
J/i P, are unchanged. In particular, in the special case when V = /C is a field, by 

Proposition 2.12, the spaces J-{H),'H{H), and the relation J/i •e-^ P, are independent 
of the minimal fractional decomposition of H. 

In the local case, when H e Mat^x^ V[(5] is a (local) Hamiltonian structure on V, then 
J/ e J^{H) = V/dV and P e -H{H) = H{d)(^lm-£^ cz are associated if and only if 

Lemma 7.3. (a) The space is a Lie algebra with bracket (6.12), and T-L{H) cz is 
its subalgebra. 

(b) We have a representation (p of the Lie algebra on the space V/dV given by 

*(f)(S/.)=|f-|. 

and the subspace J'{H) cz V/dV is preserved by the action of the Lie subalgebra 
n{H) cz 

(c) If \h -(-^ P and \h -^-^ Q for some \h e F{H), then the action of P,Q e 'H{H)on 
J~{H) is the same: 
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Proof. It follows immediately from [BDSK09, Lem.4.7-8]. □ 

Thanks to Lemma 7.3, we have a well-defined map {• , : J'{H) x J-{H) — > J-{H) 
given by 

(7.1) (5/.b)« = Jp | [-\l.M,nl), 

where P e n{H) is such that J/ ^ P- 

Proposition 7.4. (a) The bracket (7.1) is a Lie algebra bracket on the space of Haniil- 
tonian functionals J-{H) . 

(b) The Lie algebra action ofl-L{H) on J-{H) is by derivations of the Lie bracket (7.1). 

(c) The subspace 

A{H) = {(j/,p) e :f{h) X n{H) 1 5/ p] 

is a subalgebra of the direct product of Lie algebras J~{Li) x T-i{H). 

Proof. It follows immediately from [BDSK09, Prop. 4. 9, Rem. 4. 6]. □ 

A Hamiltonian equation associated to the Hamiltonian structure H and to the Hamil- 
tonian functional J/i e J^{H), with an associated Hamiltonian vector field P e 1-L{L[), is, 
by definition, the following evolution equation on the variables u = {ui)^^j- 



By the chain rule, any element / e V evolves according to the equation 

dt ^r^^ ''f),s^y 



and, integrating by parts, a local functional J/ e V/dV evolves according to 

1^ 



4/^ fp.^ 

dt J Su 



An integral of motion for the Hamiltonian equation (7.2) is a Hamiltonian functional 
^/ e T{H) such that 

i.e. ^/ lies in the centralizer of J/i in the Lie algebra T[H). The Hamiltonian equation 
(7.2) is said to be integrable if there is an infinite sequence of pairs (^/in,Pn) e ^{H) x 

T-L{H), n ^ 0, such that J/io = S^i ^Po = we have J/i^ — * Pn for every n e Z+, the 
sequence {\hn}ne'L+ spans an infinite-dimensional abelian subalgebra of the Lie algebra 
J- (Li), and the sequence {Pn}neZ+ spans an infinite-dimensional abelian subalgebra of 
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the Lie algebra 7i{H). (Equivalently, if there exists an abelian subalgebra f) of the Lie 
algebra A{H) defined in Proposition 7.4(c), such that both its canonical projections 7ri(f)) 
and 7r2(f)) in J-{H) and T-i{H) are infinite dimensional, and J/i e 7ri(f)).) In this case, we 
have an integrahle hierarchy of Hamiltonian equations 



du 
dtn 



Pn, neZ+. 



7.2 The Lenard-Magri scheme of integrabihty 

Throughout the section, we let H,K e Mat^x£V(c') be a pair of compatible Hamilto- 
nian structures on the algebra of differential functions V. We also let H = AB^^ and 
K = CD-^, with A,B,C,D e Mat£x£V[5] and detB + 0,deti:' 7^ 0, be their minimal 
fractional decomposition, and let Ca^b, C,c,D ^ V®^©V^ be the corresponding compatible 
Dirac structures (see Theorems 6.11 and 6.16). 

The Lenard-Magri scheme of integrahility consists in finding sequences 

= S/io, \hi, \h2, • • • e HH) n F{K) , 



Pn, Pi , Po 



en{H) nn{K). 



such that 
(7.3) 



P 



n-l 





Explicitly, diagram (7.3) holds if and only if there exists a sequence 
such that the following equations hold: 



(7.4) B{d)F^i 
for all n ^ 0. 



, Cid)F2n = A{d)F2n^l = Pn , B{d)F2n+l = D{d)F, 



2n 



5hn_ 

6u 



Proposition 7.5. (a) If the sequence of pairs (5^n,-Pn); ^ n ^ + 1, where S^hn e 



F{H) n F{K), Pn e n{H) n H{K), satisfies \hn-i 
then 



H 



Pn iK forQ^n^N, 



(7.5) 

(h) We have 
(7.6) 



[Prr„ Pn] e Ker 5* n Ker D* for all m,n ^ N . 
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Proof. For m = n equation (7.5) holds trivially, by the skew symmetry of the Lie brackets 
{• , and {• , ■}k- Assuming m < n, we prove (7.5) by induction on n — m. We have, 
by the assumption on the ^h^s and P„'s and the definition (7.1) of the Lie bracket on 

C Sh 

{^hrrij^hn} = Prn+1 ' = {S^rn+l, S^n} = 0, 

by the inductive assumption. Similarly, 

By Theorem 6.11, Ca,b and Cc,d are Dirac structures, hence they are closed under the 
Courant-Dorfman product (6.11). By formula (6.21) and the assumption on the ^h^s 

and PnS, we have that © P„ e £^a,b, and © P„ e Cc,d, for every ^ n ^ N . 

It then follows, by equation (6.13) and formulas (6.12) and (7.1), that 



Hence, by (7.5), we get that 0®[Pm, Pn] £ ^a,b'^^c,d- Namely, there exist F e KerS cz 
and G e KerD c such that [Pm,^n] = AF = CG e A{KeT B) n C(KerL>). To 
conclude, we finally observe that, by skewadjointness of H and K, we have B*A = —A*B 
and D*C = —C*D, which immediately implies AiKex B) a Kei B* and C(KerD) c 
KerD*. □ 

Corollary 7.6. Suppose that the sequences {J/in}nGZ+i {Pn]ne'L+, span infinite dimen- 
sional subspaces ofV/dV and respectively, and satisfy conditions (7.3) (where ^h-i = 
0). Suppose, moreover, that Ker B* nKer D* = 0. Then, the hierarchy of bi-Hamiltonian 
equations 

— = Pn , neZ+. 
dtn 

is integrable. 

Proof. It is immediate from Proposition 7.5 and the definition of integrability. □ 

Next, we want to discuss in which situations one can apply successfully the Lenard- 
Magri scheme. In order to do so, recall that an algebra of differential functions V is called 
normal [BDSK09] if -A^iVm^i) = Vm,i for every m G i ^ I, where 

du- ' ' ' 



VmA = \feV 



= for (n,j) > {m,i) (in lexicographic order) [ . 



For example, the algebra of differential polynomials in i variables is normal, and any 
algebra of differential functions V can be included in a normal one. 
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Theorem 7.7. LetV be a normal algebra of differential functions. Let H , K eMa.texeV{d) 
be a pair of compatible non-local Hamiltonian structures on V and assume that det K ^ 

0. Let H = AB^^ , K = CD^^, he their minimal fractional decompositions, with 
A,B,C,De MaUxiV[d]. Let N ^ 0, and letO = J/iq, J/ii, • • • , J/iTV e F[H)nF{K) 

and Po,Pi,...,Pn e n{H) n ?^(A') be such that ^h^-i ^ Pn ^ for Q ^ n ^ N . 
Suppose, moreover, that the following orthogonallity conditions hold: 

(Spanc{^}^)lo)^ ^I^C- , (Spanc{i^4^^o)^ cImS, 

where the orthogonal complement is with respect to the pairing V®^ x ^ V/dV, defined 
by {F, P) ^ S^F ■ P. Then, there exist infinite sequences {^hn}nez+ ^[H] n F{K) and 
{-Pn}nsZ+ T~^{H) n HlK), extending the given finite sequences, satisfying conditions 
(7.3) . In other words, the Lenard-Magri scheme of integrability can be applied. 

Proof. By Proposition 7.5, we have 

l^-^-Pn = {lhN,<Sh^}K = 0, 

for all ?i = 0, . . . , N. Hence, by the second orthogonality condition we get that = 
B{d)X for some X e V®^. Let Pn+i = A{d)X e V^. Again by Proposition 7.5, we have 

^~Iu ' ^^^^ ^ ^^^^^ ^ ^ ' 

for all n = 0, . . . , A''. Hence, by the first orthogonality condition we get that Pn+i = 
C{d)Y for some Y e V®^ Let Fjv+i = D{d)Y e V®^ It follows by Theorem 6.19(c) (with 
in place of Fq, in place of Fi, and F/v+i in place of F2) that F^r+i is exact, 

1. e. Fat+i = for some J/iAf+i £ V/dV. By construction, J/ijv <— Pn+i *— S^Af+i- 
Hence, we extended the sequences {\hn}n=-i and {Pn}n=o by one step. The claim follows 
by induction. □ 

Remark 7.8. Let V be an algebra of differential function in Ui, i e /. Assume that V is a 
domain, and let /C be its field of fractions. Let H = AB^^, K = CD^^ e Mat^x^ V(5) be 
skewadjoint rational matrix differential operators with coefficients in V, with A,B,C,D e 
Matixi V[d], det B, det D 0. Let = 0, Co, • • • Itv e V®^ and Po,...,Pn eV^ be such 
that 

(7.7) Cn-l = B{d)Fn , Pn = A{d)Fn , Cn = D{d)Gn , Pn = C{d)Gn , 

for ^ n ^ A^ and some Fq, . . . , F/v, Go, . . . , Gn £ V^, and 

( Span c }lo) ^ c Im C , ( Span c{Pn}Lo)^ ^ ■ 

By the same arguments as in the proof of Theorem 7.7 we can extend the given seguences 
to infinite sequences {£,n}nez+ c V®^ and {Pn}neZ+ ^ satisying conditions (7.7) for 
every n e Z+. Next, suppose that H = AB^^ and K = CD^^ are minimal fractional 
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decompositions (cf. Definition 2.11) and that H and K are compatible Hamiltonian 
structures on V. Suppose, moreover, that detK j^O and Co = §f for some J/iQ e V/dV. 
Then it follows by Theorem 6.19 that all C„,'s are closed in /C®^, i.e. D^^{d) is self-adjoint 
for every n e Z+. Hence, all the C„'s are closed in V®^. Finally, taking, if necesary, a 
normal extension V of V, one proves that all the ^^'s are exact, i.e. there exist elements 

iK e V/dV such that C„ = ^Jt' ^ e Z+, [BDSK09, Thm.3.2]. Hence, the Lenard-Magri 
scheme can be applied. 

Remark 7.9. Recall from Theorem 5.1 that, given two compatible non-local Hamiltonian 
structures H and with det(A') 7^ , we get an infinite family of compatible non-local 
Hamiltonian structures given by ijM = K, and ifW = {H o K'^f-^ o H, for s ^ 1. 
Suppose that the infinite sequences {\hn]nez+ c F{H) n T{K), {Pn}neZ+ ^ T~(-{H) n 
7i{K) satisfy the Lenard-Magri recursive relations J/in-i ^ -fn "—^ for every n ^ 

(we let 5^-1 = Jo)- Then, for every n, s e Z+ we have the relations ^hn ^ Pn+s- This 
will be proved in [DSK12]. Consequently, all the evolutionary equations ^ = Pn are 
Hamiltonian with respect to all the non-local Hamiltonian structures H^'^\ s e Z+. 

7.3 Example: NLS integrable hierarchy 

Consider the algebra of differential polynoamials in two variables R2 = F[m, v, u' , v' , . . .], 
and its extension Vu = ¥[u-^ ,v,u' ,v' , . . .]. Let also IC be the field of fractions of R2 
(which is the same as that of Vu). Consider the following pair of compatible non-local 
Hamiltonian structures with coefficients in R2 from Example 4.10: 



vd ^ o V —vd ^ o u \ / — 1 

^-1 ] +Cdl , K = { 

-uo ovud ouj ylU 

where c e F is a fixed constant. The operators H and K admit the following fractional 
decompositions: H = AB^^, K = CD^^, where C = K, D = 1, and 

^ . cd o u —v?v \ B ( ^ ^ 

cd ov + cdo {ud + 2u') y ' \ v ud +2u' 



This fractional decomposition of H is minimal over /C by Proposition 2.12 since Ker i? 
K? is spanned by ( 2 ) ' which does not lie in Ker A. 



Let \h-i = 0, 5/10 = \{u^ + v"^) e VjdVu and Pq = ^ ""^"^ 

= 5^-1 <— » Pq 5^0- Indeed, letting F_i = ^ ^'^2 ^ ^i^d ^ 7; ^ ' have 
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Note that, even though = S^hi,S^ho e R2/dR2 and Pq e R2, they are NOT associated 
withm the algebra R2, since F-i does not he in R2. 

Next, we prove the orthogonahty conditions: (-1^)''^ Im I, {Po}^ Im.{B : 

Va. The first condition . obvious. For the second. ^ ^ ( ^ ) e ts otthogon^ to P„ 

if and only if 

—vf + ug = dh e dVu ■ 

In this case, F = B{d) j^j^^ ) e ImS. 

By the observations in Remark 7.8, there exist sequences {■^n}nez+ and {Pn}neZ+ in 
satisfying condition (7.7) for all n e Z+, and all the ^„,'s are closed, in the sense that 
their Frechet derivative D^^{d) is self-adjoint for every n. 

We argue that, in fact, all these elements P„, n e Z_|_, lie in i?2- Indeed, it is not 
hard to show, using the explicit form of the matrices A and B, that, if ^n-i lies in i?2, 
then necessarily both P„ and he in i?2- Alternatively, we can exchange the roles of u 
and V to prove that all ^„'s and Pn's lie in V^, hence in i?2. First, we observe that, since 
Ker B is one dimensional over F and C = K is invertible, the flag of subspaces 



ON Uq = Spanip{^o} , Ui = Spanip{^o, Ci} , • • • , 

^ ^ Vo = Spanj.{Po} , Vi = Spanp{Po, A } , • • • c , 



associated to the sequences {Cn}neZ+ and {Pn}neZ+ satisfying condition (7.7) for every 
n e is uniquely defined. Consider now the new fractional decomposition H = AiB^^, 
where 

^ f + cd o {vd + 2v') cd ou \ ^ f vd + 2v' u 
^ y —uv^ cd o V J ^ ^ \ V 

By the same argument as above, we get that there exist sequences {Cn}neZ+ and {P^}neZ+ 
in V^, satysfying condition (7.7) with respect to the new fractional decomposition, i.e. 

en-i = Bi{d)F}, , Pi = Ai(5)F„i, en = D{d)Gi, P^ = C{d)Gl, 

for some elements F^, G\ e V^. Note that the two fractional decompositions H = AB^^ 
and H = AiB^^ are equivalent, in the sense that Ai and Bi are obtained from A and B 
by multiplication on the right by the matrix 

\ —u ^v'^ 

which is invertible in Mat£x£^[(^]- Hence, by the arguments in Remark 7.8 we get that 
both the sequences {^n}nEZ+, {Pn}neZ+ c V^, and the sequences {C^}nGZ+, {Pn}neZ+ c 
V^,, considered as sequences in /C^, satisfy condition (7.7) for every n e Z+. Therefore, by 
the uniqueness of the flag (7.8), we get that the subspaces C/„ = Spaupj^O) ■ • • iCn} 
and = Spauifl^g, . . . ,^^} cz must be the same subspace of /C^, and similarly the 
subspaces Vn = Span^{PQ, . . . , P„} cz and V^l = Spanip{PQ^, . . . , P^} cz must be the 
same subspace of /C^. In conclusion, all ^„'s and P^s lie in n = i?2- 
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Since R2 is a normal algebra of differential functions, it follows from [BDSK09, 
Thm.4.13] that all ^ n's are exact, i.e. there exists elements J/in £ V/dV, n £ Z-f, such 

that ^ri = "IIt"- Hence, the Lenard scheme can be applied. 

It is not hard to show, by induction on n, that both entries of e i?| have differential 
order n. Therefore, the vectors are linearly independent. It follows that the elements 
S^hn e V/5V, n e Z+ and the elements Pn = K^n ^ R^-, n e Z+, are linearly independent 
as well. Since D = I, we can apply Proposiion 7.5 and Corollary 7.6 to deduce that we 
have an infinite hierarchy of compatible integrable equations ^ = P„, n e Z+, for which 
S^hn e J~{H) n J-{K) are integrals of motion. 

We can compute the first few equations of this integrable hierarchy, called the non- 
linear Schroedinger (NLS) hierarchy: 

d f u \ _ f —V 
dtQ \ V J \ u 

dti\v ) ~ ^\ v' 
d ( u \ _ ( c^v" + f f (n^ + v^) 



dt2 \ V J \ -c^u" - ^u{v? + 
d f u\ ( -c^u'" - \c^{u^ + v'^)u' 



dt3 \ V J \ -c^v'" - \c^{u^ + v'^)v' 

The third equation is called the (coupled) NLS equation. The first four integrals of 
motion are 



= J _ c^uv'" - "^{u'v' - v\')) . 



Of course this is a very well known hierarchy, studied by many different methods, see 
e.g. [TF86, Dor93, BDSK09]. The approach of the last two papers, based on a Dirac 
structure, is close to ours (even though in [Dor93] the proposed Dirac structure is not 
quite a Dirac structure). 
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